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ABSTRACT 

We report the discovery of a double-lined, spectroscopic, eclipsing binary in the 
Orion star- forming region. We analyze the system spectroscopically and photomet- 
rically to empirically determine precise, distance-independent masses, radii, effective 
temperatures, and luminosities for both components. The measured masses for the pri- 
mary and secondary, accurate to ~ 1%, are 1.01 M and 0.73 M Q , respectively; thus the 
primary is a definitive pre-main-sequence solar analog, and the secondary is the lowest- 
mass star yet discovered among pre-main-sequence eclipsing binary systems. We use 
these fundamental measurements to test the predictions of pre-main-sequence stellar 
evolutionary tracks. None of the models we examined correctly predict the masses of the 
two components simultaneously, and we implicate differences between the theoretical 
and empirical effective temperature scales for this failing. All of the models predict the 
observed slope of the mass-radius relationship reasonably well, though the observations 
tend to favor models with low convection efficiencies. Indeed, considering our newly 
determined mass measurements together with other dynamical mass measurements of 
pre-main-sequence stars in the literature, as well as measurements of Li abundances 
in these stars, we show that the data strongly favor evolutionary models with ineffi- 
cient convection in the stellar interior, even though such models cannot reproduce the 
properties of the present-day Sun. 
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1. Introduction 

Pre-main-sequence stellar evolutionary models are fundamental to our paradigm of star for- 
mation and early stellar evolution. These models are the means by which masses and ages are 
assigned to young stars and brown dwarfs, and are therefore central to our understanding of phys- 
ical processes that depend, directly or indirectly, on knowledge of stellar masses and ages. Indeed, 
pre-main-sequence (PMS) stellar evolutionary models touch upon the most basic questions in star 
formation research: the nature and origin of the initial mass function; the timescale for circumstellar 
disk evolution and planet formation; and the initial distribution of stellar angular momentum, its 
evolution in time, and the dependence of that evolution on stellar mass, accretion, and other mass- 
dependent properties. Thus, the absolute calibration of PMS evolutionary models with respect to 
stellar mass is of foundational importance to star formation research. 

Empirical mass determinations exist in the literature for 15 individual PMS stars (Popper 
1980; Andersen 1991; Casey et al. 1998; Simon, Dutrey, & Guilloteau 2000; Covino et al. 2000; 
Steffen et al. 2001; Alencar et al. 2003), with only 6 of these in the mass domain of 1 M or below. 
Moreover, many of these measurements, while of high precision, potentially suffer from systematic 
errors (primarily distance uncertainties) that limit their accuracy. Indeed, there is no PMS star 
with M < 1 M whose mass has been measured with an accuracy of better than a few percent. 

Consequently, the absolute mass calibration of theoretical PMS evolutionary tracks remains 
weakly constrained by observations. The dearth of accurate measurements below 1 M Q is partic- 
ularly salient because low-mass stars dominate studies of young stellar populations, and it is in 
this mass domain that key, outstanding questions remain with respect to, e.g., the nature of the 
initial mass function (Hillenbrand 1997; Muench et al. 2002) and the evolution of stellar angular 
momentum (SMMV; Herbst et al. 2002; Rebull 2001). 

Eclipsing binary stars provide a uniquely powerful means of determining absolute stellar masses 
because the mass determinations are distance-independent. Eclipsing binaries also yield direct 
measurement of the stellar radii, permitting the determination of stellar luminosities that are also 
distance independent. To date, only four eclipsing binaries with PMS components have been 
reported in the literature (Andersen 1991; Popper 1980; Casey et al. 1998; Covino et al. 2000). 

Here we report the discovery of a previously unknown, double-lined, spectroscopic, eclipsing 
binary system in the Orion star-forming region, V1174 Ori. We combine time-series photometry 
with high-resolution spectroscopy to determine all of the fundamental physical properties of the 
system and of its stellar components. We then use these measurements to test the predictions of a 
variety of PMS stellar evolutionary models. 

We begin by describing the discovery of V1174 Ori and its identification as a PMS eclipsing 
binary in §2. In §3, we present our photometric and spectroscopic observations and the methods 
used in their reduction. We analyze these observations in §4 to determine the system ephemeris 
and to derive a precise double-lined orbit solution, providing absolute stellar masses. We conduct 
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a thorough spectroscopic analysis that includes determination of the stellar rotational velocities, 
measurement of the primary-to-secondary flux ratio, spectral classification and determination of 
effective temperatures, measurement of Li abundances, and measurement of the equivalent width 
of Ha emission. We also present a full analysis of our multi-band time-series photometry using 
synthetic light-curve fitting to determine precise, absolute stellar dimensions. The results of our 
analysis are presented in §5, including an independent determination of the distance to VI 174 Ori. 

In §6 we use our empirically determined stellar masses, radii, effective temperatures, and lumi- 
nosities to test a variety of PMS stellar evolutionary tracks possessing a variety of characteristics, 
including different atmosphere models and treatment of convection. We assess the performance of 
the theoretical models both against the components of VI 174 Ori alone as well as against the en- 
semble of all PMS stars with empirical mass determinations in the literature. We consider how the 
tracks compare quantitatively to the data in both high- and low-mass regimes, in both binary and 
single systems, and in both the H-R diagram and mass-radius observational planes. Moreover, we 
examine the general behavior of the tracks to suggest, in qualitative terms, how their performance 
against the data might be improved. 

In addition, we assay the ensemble of empirical mass determinations independent of the models 
to explore the degree of internal consistency in the data themselves, and we show how Li abundances 
can be used in some cases to resolve apparent contradictions. Finally, we discuss more general 
insights that can be gleaned from the Li abundances, arguing that the data strongly favor models 
with inefficient convection in the stellar interior during the PMS stage of evolution, even though 
such models cannot reproduce the properties of the present-day Sun. We present a summary of 
our conclusions in §7. 



2. Discovery 

Here we describe the discovery of VI 174 Ori as a PMS eclipsing binary system. Its position is 
a = 05:34:27.85, 6 = -05:41:37.8 (J2000.0). 



2.1. Data from the Literature 

VI 174 Ori has been included in a number of previous surveys of stars in the vicinity of the 
Orion Nebula Cluster (ONC), dating back to the catalog of Parenago (1954) (star # 1478). The 
General Catalog of Variable Stars, whose designation we use in this paper, reports a range of 1.1 
mag (in the photographic system) between its high and low states. 

More recently, the ROSAT observatory included V1174 Ori (designation ROS-ORI 166) in 
a deep PSPC pointing centered on the Trapezium region. The analysis of these data by Geier, 
Wendker, & Wisotzki (1995) assigns a net of 75 X-ray photons to V1174 Ori, for a count rate 
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of 17.02 cts/ks. This corresponds to an X-ray luminosity of Lx ~ 1 x 10 30 erg/s in the ROSAT 
bandpass (0.07 to 2.4 keV), assuming the ONC distance to V1174 Ori of 470 pc (e.g. Genzel et al. 
(1981)). V1174 Ori has not been included in any of the Chandra observations of the ONC currently 
in the public archive, but this X-ray luminosity is fairly typical among the PMS stars that have 
been observed by Chandra (e.g. Feigelson et al. (2002)). 

V1174 Ori was observed in the K-band survey of the Orion A cloud by Ali & Depoy (1995) (# 
2176 in their catalog), but their reported -fT-band magnitude of 9.00 is at the bright limit of their 
survey, making this measurement suspect. More recently, Carpenter, Hillenbrand, & Skrutskie 
(2001) presented JHK light curves of V1174 Ori (star number 4008 in their catalog), as part of 
their analysis of time-series photometry of the Orion star-forming region, conducted as an extension 
of the 2MASS project 7 . Their light curves show rms scatters of 0.08-0.09 mag. In addition to low- 
level variability, their light curves display evidence for a possible eclipse event in the form of a single 
observation (seen in all three filters) that is 0.3-0.4 mag dimmer than the mean brightness level. 

Finally, VI 174 Ori was included in the study of Rebull et al. (2000) (# 1-665 in their catalog), 
whose UVI photometry places V1174 Ori on the locus of Orion PMS stars in the (U — V) vs. 
(V — I) color-magnitude diagram. They report a spectral type for VI 174 Ori of M0 on the basis 
of low-resolution spectroscopy, but derive extinctions of Ay = —0.50 and Aj = —0.31, indicating 
potential problems with the spectral classification or photometry. While V1174 Ori was included 
in the region covered by the time-series study of Rebull (2001), they do not report it as a periodic 
variable. VI 174 Ori was in eclipse during at least 7 of the 30 nights covered by their observations 
(assuming the ephemeris we adopt in §4) , but with only one observation per night their light curves 
were presumably too sparse to detect periodic behavior (if they detected the eclipses at all). 

2.2. Identification as a PMS Double-Lined Eclipsing Binary 

We first identified V1174 Ori as a candidate eclipsing binary from an /-band light curve 
obtained by us in 1994 December as part of a photometric study of stellar rotation in the ONC 
(Stassun et al. 1999; hereinafter SMMV). The discovery light curve was by itself very compelling, 
showing eclipse-like variability with a period of s=s 2.63 days. 

V1174 Ori was also included in several pointings during our spectroscopic survey of the region 
with the WIYN 8 multi-object spectrograph, which showed V1174 Ori to have strong Li absorption 
(an indicator of stellar youth) and revealed it to be a single-lined spectroscopic binary. We confirmed 
the spectroscopic presence of the secondary in a spectrum kindly obtained by G. Basri with the 
Keck HIRES spectrograph, which also clearly revealed the presence of Li absorption due to the 



7 See http://www.ipac.caltech.edu/2mass for information about 2MASS 

8 The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, Yale 
University, and the National Optical Astronomy Observatories. 
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secondary (Fig. 1). This spectrum also allowed us to determine that the primary was of roughly 
mid-K spectral type. 

3. Observations and Data Reduction 

Having identified V1174 Ori as a PMS eclipsing binary and double-lined spectroscopic bi- 
nary, we obtained multiple high-S/N, high-resolution spectra with which to determine a precise, 
double-lined orbit solution. We also obtained high-precision, high-cadence light curves at multiple 
wavelengths in order to (a) firmly establish the system ephemeris, and to (b) determine, via light- 
curve synthesis and modeling, such parameters as the orbital inclination, component radii, and 
ratio of component effective temperatures. The one remaining piece of information needed to fully 
establish all system parameters is the effective temperature of one of the components, a datum also 
provided by our high-resolution spectroscopy. 

Here we present the photometric and spectroscopic data we obtained, and describe the proce- 
dures we employed in their reduction. 

3.1. Photometry 

We observed V1174 Ori photometrically with CCDs at lm-class telescopes at 6 separate epochs 
spaced over 8 years. Table 1 summarizes the time-series photometry we obtained. The individual 
measurements are available electronically and are shown in Figs. 2, 4, and 6. 

All CCD data were processed in the manner described in SMMV. Briefly, we reduced our 
CCD frames using the IRAF 9 CCDRED or XCCDRED (for CTIO quad-readout data) packages, and 
performed aperture photometry on all stars in the field of view using the IRAF APPHOT package. 
We then applied an algorithm similar to that of Honeycutt (1992) to our raw light curves to remove 
non-cosmic frame-to- frame photometric variations (due to changes in, e.g., seeing, sky brightness, 
atmospheric transparency). In all cases, the resulting photometric precision in the light curves of 
V1174 Ori is limited by systematics (flat-fielding, PSF variations, etc.) at the ~ 0.01 mag level, 
which we determine from inspection of non-variable stars in the field of comparable brightness. 

In addition to these differential time-series observations, we obtained calibrated BVI photom- 
etry of VI 174 Ori on 1995 Feb 02 at the USNO lm telescope. These data were obtained out of 
eclipse, at orbital phase 0.35, based on the ephemeris presented in §4.1. Calibrated UVI pho- 
tometry of V1174 Ori has been reported by Rebull et al. (2000), and Carpenter, Hillenbrand, Sz 
Skrutskie (2001) report JHK measurements. These measurements are summarized in Table 2. 



9 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 
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Rebull et al. (2000) do not report the precise times of their observations. However, the V — I 
color reported by them is identical to that measured by us, which leads us to believe that their 
observations were also obtained out of eclipse. Carpenter, Hillenbrand, & Skrutskie (2001) report 
mean JHK magnitudes and colors, averaged over all of their observations which span 36 days 
(excluding the two non-contiguous measurements obtained in 1998 Mar and 2000 Feb). The J 
magnitude we adopt is the mean of their non-eclipse observations. The J — H and H — K colors 
do not show significant variations in their observations and so we do not adjust the values reported 
by those authors. 



3.2. Spectroscopy 

In order to derive a precise, double-lined orbit solution for V1174 Ori, we obtained 15 high 
S/N observations in queue observing mode with the High Resolution Spectrograph (HRS) on the 
Hobby Eberly Telescope 10 (HET) on 15 different nights between 2001 Nov 22 and 2002 Feb 23. 
The spectrograph setup we used yielded a resolving power of R k, 30000, with wavelength coverage 
from 5095A to 8860A, and centered at 6948A. The 51 echelle orders in each spectrum are imaged 
onto two CCDs, with the "blue" chip imaging the 32 orders from 5095A to 6803A and the "red" 
chip imaging the remainder. The spectrograph was also set up with a single sky fiber so that 
simultaneous sky spectra could be obtained with each observation. The CCD images were binned 
on readout such that a spectral resolution element corresponds to approximately 3 pixels in the 
extracted spectra. 

The exposure times were initially 3 x 600 = 1800 sec, and later increased to 3 x 800 = 2400 sec, 
and we combined the individual exposures in each set with cosmic-ray rejection. Each observation 
of V1174 Ori was bracketed by an observation of a ThAr lamp. In addition, all but one observation 
was accompanied by an observation of a bright radial- velocity standard star (typically obtained 
shortly after the VI 174 Ori exposure) in order to closely monitor and correct for instrumental 
drifts in the radial velocities. 

In addition, we obtained a single, very high S/N (~ 250) observation of a late-type radial- 
velocity standard star for use as a radial- velocity template in our cross-correlation analysis (§4.2.2). 
Once we had established the approximate spectral types of the components of V1174 Ori (§4.3), 
we also obtained single observations of two appropriate spectral-type standard stars, of types K3 
and K7 11 . Table 3 summarizes our spectroscopic observations with HET HRS. 



The Hobby-Eberly Telescope is operated by McDonald Observatory on behalf of The University of Texas at 
Austin, the Pennsylvania State University, Stanford University, Ludwig-Maximilians-Universitat Miinchen, and 
Georg-August-Universitat Gottingen. The observations described here were obtained through community access 
made possible by NO AO. 

11 More careful analysis of the V1174 Ori spectra after the observations were obtained indicates that the components 
of V1174 Ori have spectral types of K4.5 and ~M1.5 (see §4.3.3), but this does not significantly affect our analysis. 
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Table 1. Photometric time-series observations of V1174 Ori 



UT Dates 


N (nights) 


Range of Julian Dates 
(2400000+) 


Observatory 


Filter 


# of 
obs. 


Ref. 


1994 Dec 11-27 


16 


49698.340-49714.491 


KPNO 0.9m, USNO lm, Wise lm 


Ic 


128 


SMMV 


2001 Jan 20-27 


8 


51929.588-51936.775 


CTIO 0.9m 


Ic 


152 


this study 






51929.608-51936.768 




V 


144 


this study 






51929.590-51936.776 




Str. y 


158 


this study 






51929.594-51936.783 




Str. b 


154 


this study 






51929.623-51936.735 




Str. v 


130 


this study 


2001 Jan 31 


1 


51940.591-51940.819 


USNO lm 


Ic 


20 


this study 






51940.589-51940.817 




V 


20 


this study 


2001 Nov 14-24 


9 


52227.746-52237.996 


WIYN 0.9m 


Ic 


78 


this study 


2002 Nov 17-26 


10 


52595.745-52605.000 




Ic 


110 


this study 






52595.802-52605.004 




V 


96 


this study 






52595.806-52604.984 




B 


96 


this study 


2002 Dec 05-16 


11 


52613.810-52624.946 




Ic 


24 


this study 



Table 2. Calibrated photometry of VI 174 Ori 



Filtcr(s) 


Mag./Color 
observed 


Mag./Color 
calculated 81 


UT Date(s) 


Ref. 


V 


13.95 ±0.03 


13.95 


1995 Feb 02 


1 


B-V 


1.25 ±0.02 


1.24 




1 


V -I 


1.62 ±0.02 


1.48 




1 




1.61 ±0.02 




1998 Jan 


2 


U - V 


2.22 ±0.02 


2.33 




2 


J 1 


11.22 ±0.05 


11.52 


2000 Mar 04 - Apr 08 2 


3 


J - H 


0.65 ±0.02 


0.65 




3 


H -K 


0.19 ±0.03 


0.14 




3 



References. - 1: this study, 2: Rebull et al. (2000), 3: Carpenter, 
Hillenbrand, & Skrutskic (2001) 

a Assuming main-sequence colors, bolometric corrections, and Ry = 3.12. 
See §5.1. 

1 Value adjusted from that reported in Carpenter, Hillenbrand, & Skrut- 
skic (2001), in order to account for eclipses. 

2 Contiguous portion of light curve only. 
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Table 3. Log of HET HRS Observations 



II 

# 


U 1 Date 


T T TT~vl 


Object 


Exp. Time (s) 


o/iN 


Comments 




o 


2001 Oct 29 


2212.74347 


HD 18884 


30 


245 


R.V. template, SpT = 


M2 III 


1 


2001 Nov 22 


2235.85217 


V1174 Ori 


1800 


85 










2235.87681 


HD 26162 


5 


120 


R.V. standard, SpT = 


K2 III 


2 


2001 Nov 23 


2236.84827 


V1174 Ori 


1800 


80 










2236.89321 


HD 26162 


15 


145 






3 


2001 Dec 10 


2253.80559 


V1174 Ori 


1800 


60 










2253.82748 


HD 26162 


10 


85 






4 


2001 Dec 18 


2261.78887 


VI 174 Ori 


1800 


40 










2261.80695 


HD 26162 


10 


90 






5 


2001 Dec 21 


2264.77976 


VI 174 Ori 


1800 


80 










2264.79712 


HD 26162 


10 


95 






6 


2001 Dec 22 


2265.76822 


V1174 Ori 


1800 


75 










2265.83415 


HD 26162 


10 


95 






7 


2001 Dec 23 


2266.76687 


V1174 Ori 


1800 


105 










2266.81776 


HD 26162 


10 


120 






8 


2002 Jan 26 


2300.67336 


V1174 Ori 


2400 


125 






9 


2002 Jan 27 


2301.66941 


V1174 Ori 


2400 


95 










2301.72972 


HD 26162 


10 


100 






10 


2002 Jan 28 


2302.66839 


V1174 Ori 


2400 


105 










2302.72798 


HD 26162 


10 


140 






11 


2002 Feb 07 


2312.64124 


V1174 Ori 


2400 


100 










2312.70694 


HD 26162 


10 


105 






12 


2002 Feb 09 


2314.63298 


V1174 Ori 


2400 


120 










2314.69301 


HD 26162 


10 


160 






13 


2002 Feb 13 


2318.63282 


V1174 Ori 


2400 


95 










2318.68518 


HD 26162 


10 


130 






14 


2002 Feb 21 


2326.61211 


V1174 Ori 


2400 


90 










2326.62791 


HD 26162 


10 


115 






15 


2002 Feb 23 


2328.60823 


V1174 Ori 


1800 


100 










2328.62179 


HD 26162 


10 


110 






16 


2002 Mar 07 


2340.93213 


HD 110463 


180 


140 


SpT standard, SpT = 


K3 V 






2340.64503 


HD 237903 


240 


160 


SpT standard, SpT = 


K7 V 



heliocentric Julian Date (2450000+) 
Signal-to-noise is per 3-pixel resolution element, measured near 6500A. 
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All spectra were processed with the IRAF CCDPRDC and ECHELLE packages. The steps involved 
in producing the final spectra include bias subtraction, flat-fielding, bad-column interpolation, 
order tracing (for both object and sky apertures), sky subtraction, wavelength calibration, and he- 
liocentric velocity correction. Extraction of the spectra from the raw data was not straightforward, 
particularly with respect to order tracing which required a high degree of manual intervention. The 
orders are highly angled on the CCD, have flat-topped profiles in the cross-dispersion direction, and 
are separated from the sky orders by only a few pixels. We thus encountered difficulties with the 
automated order-finding and tracing frequently "skipping" between the object and sky apertures. 
However, we found that these problems could be circumvented manually. 

ThAr arc spectra were extracted similarly. We found that relatively high-order polynomials 
both along and across the orders were required to provide a good mapping between pixel, order 
number, and wavelength. In the end we adopted third-order and seventh-order polynomial fits along 
and across the orders, respectively. The wavelength solutions were applied to the object spectra by 
taking an average of the ThAr spectra taken before and after the object spectra, weighted by their 
proximity in time to the object spectra. 

Heliocentric velocity and date corrections were computed via the IRAF RV package. We care- 
fully assigned mid-exposure observation times to each combined spectrum, taking into account the 
flux levels of the 3 sub-exposures (which occasionally differed from one another) and their time 
spacing relative to one another. 

As ours was the first community program to use HET HRS, we unfortunately encountered an 
unforseen problem that rendered the "red" orders in our spectra unusable. The "red" HRS CCD in 
use at the time of our observations suffers from severe fringing, which would normally be removed 
by a spectrum of a bright continuum source such as an internal lamp. However, in the case of 
the HRS, the fringing pattern produced by the internal flat-field calibration lamp is different from 
the pattern present in the data. This is in part due to the fact that a different fiber is used for 
calibration purposes, which results in the flat-field beam having a different //-ratio as compared 
to the science-object beam (M. Shetrone, priv. comm.). Thus in our spectroscopic analyses we use 
only the orders on the "blue" CCD, which show no such fringing problems. 

4. Analysis 

4.1. Eclipse timings: System ephemeris 

We use eclipse timing measurements from our photometric light curves (§3.1) to determine a 
precise period and ephemeris for VI 174 Ori. Our /-band data, which span a total of 2927 days, 
include 7 primary and 5 secondary eclipses. Some of these eclipses were also observed in other filters, 
providing additional (though not fully independent) eclipse timing information. We measure the 
time of each eclipse minimum by fitting a gaussian to the points in the light curve in and around 
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the eclipse. These eclipse timings and their uncertainties are reported in Table 4. 

We searched the combined /-band lightcurves for a period via the Phase-Dispersion Minimiza- 
tion (PDM) technique (Stellingwerf 1978), and found a best period of 2.61156 days. However, visual 
inspection of the combined data folded on this period revealed that the algorithm had chosen a beat 
period whereby the primary eclipses from the 2001 data were overlaid on the secondary eclipses 
from the 1994 data. Thus we refined the period by searching our radial-velocity measurements 
(§4.2.2) for periods near this value. This yielded a best period of 2.63469 ± 0.00012 days. A final 
visual inspection of the light curve folded on this period suggested a small manual refinement to 
P = 2.634727 ± 0.000004 days. 

In Fig. 8 we show the O — C phase residuals of the eclipse timings from Table 4 using this 
period. To show the residuals of the secondary eclipse together with those of the primary, we have 
subtracted 0.5 phase from the secondary residuals, which assumes a circular orbit. These residuals 
have an r.m.s. of 0.0015 in phase (5.5 min), and show no trends over the 2900-day span (some 1100 
cycles) of the data. While there appears to be a weak effect in the Jan 2001 data for the primary 
eclipse (O — C) values to be systematically negative and for the secondary eclipse (O — C) values 
to be systematically positive, it should be borne in mind that the eclipse timings from the different 
filters in Jan 2001 are not truly independent, as these light curves share a common sampling pattern. 
The eclipse timings are therefore consistent with the assumption of a circular orbit. 

We thus adopt the following ephemeris for the remainder of our analysis: 

Min I : HJD 2, 449, 703.7504(15) + 2.634727(4)^, 

where E is the epoch number, and the zero-point corresponds to primary eclipse (i.e. the deeper 
minimum). 

4.2. Radial velocities and orbit solution 

We use our HET HRS observations (Table 3) to measure precise radial velocities for both 
components of the V1174 Ori binary. Here we describe our procedure for measuring radial velocities, 
and discuss our assessment of potential systematic effects that may degrade the accuracy of the 
masses that we derive from them. 



^.2.1. Radial velocities 

We measure radial velocities by cross-correlating each VI 174 Ori spectrum with a high-S/N 
spectrum of the radial- velocity standard HD 18884 (see Table 3), for which we adopt the heliocentric 
radial velocity of —26.1 km s _1 as determined by Udry et al. (1999), whose long-term monitoring 
of this star shows it to be non-variable to a precision level of 0.3 km/s over nearly 300 observations. 
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Table 4. Times of Eclipse Minima 



UT Date 


HJD 1 


Filter 


Eclipse type 


O-C 

(phase) 


1994 


Dec 11 


4969S 4867 + 0102 


In 

i C 


Primary 


0.0022 


1994 


Dec 13 


49699.7974 ± 0.0193 


Ic 


Secondary 


-0.0004 


1994 


Dec 15 


49702.4341 ± 0.0027 


Ic 


Secondary 


0.0004 


1994 


Dec 17 


49703.7433 ± 0.0010 


Ic 


Primary 


-0.0027 


1994 


Dec 21 


49707.7129 ± 0.0158 


Ic 


Secondary 


0.0040 


1994 


Dec 27 


49714.2850 ± 0.0011 


Ic 


Primary 


-0.0016 


2001 


Jan 23 


51932.7275 ± 0.0004 


Ic 


Primary 


-0.0007 






51932.7284 ± 0.0003 


V 




-0.0004 






51932.7274 ± 0.0003 


y 




-0.0008 






51932.7276 ± 0.0002 


b 




-0.0007 






51932.7278 ± 0.0003 


V 




-0.0006 






51932.72774 ± 0.00020 


mcan a 




-0.00065 


2001 


Jan 27 


51936.6835 ± 0.0009 


Ic 


Secondary 


0.0008 






51936.6833 ± 0.0014 


V 




0.0007 






51936.6859 ± 0.0013 


y 




0.0017 






51936.6830 ± 0.0011 


b 




0.0006 






51936.6775 ± 0.0026 


V 




-0.0015 






51936.68264 ± 0.00155 


mcan a 




0.00042 


2001 


Jan 31 


51940.6324 ± 0.0004 


Ic 


Primary 


-0.0005 






51940.6331 ± 0.0003 


V 




-0.0002 






51940.63275 ± 0.00049 


mcan a 




-0.00033 


2001 


Nov 14 


52227.8162 ± 0.0011 


Ic 


Primary 


-0.0010 


2001 


Nov 18 


52231.7679 ± 0.0032 


Ic 


Secondary 


-0.0012 


2001 


Nov 22 


52235.7250 ± 0.0026 


Ic 


Primary 


0.0007 



heliocentric Julian Date (2400000+) 

a Mcan of the measurements obtained from the individual filters. 



Note. — O—C values are with respect to the ephemcris for primary minimum 
HJD 2,449,703.7504 + 2.634727S (see text). 
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We chose this late-type (M2) template in the hopes of maximizing the cross-correlation peak of the 
faint, late-type secondary star. 

We cross-correlate each of the 32 spectral orders in a given spectrum separately, and mea- 
sure the centroid of the primary and secondary star peaks (see Fig. 9a). For each spectrum we 
thus in principle have 32 separate (but not independent) radial-velocity determinations for each 
of the primary and secondary. In practice, however, some orders consistently produce poor radial 
velocities, typically due to the presence of very broad lines (e.g. the Na D lines) which result in 
cross-correlation peaks that do not lend themselves well to centroiding. In addition, we were unable 
to securely identify the peak corresponding to the secondary star in most orders, due to other low- 
level structure in the cross-correlation function. To be conservative, we thus typically only accepted 
a few very good orders, orders where the spectral features of the secondary are particularly strong. 
Finally, we found that it was necessary to avoid spectral regions with strong nebular emission, 
which often subtracted poorly. By carefully selecting the spectral regions to be used in each order 
for each spectrum, we were able to maximize the number of orders yielding radial velocities. 

For the secondary star, we use all orders that unambiguously show a secondary peak to define 
a mean radial velocity in each spectrum. With the primary star, for which most of the spectral 
orders yield a velocity measurement, we can be even more discriminating to maximize our precision. 
As shown in Fig. 9b, the scatter in the primary radial velocities produced by the different spectral 
orders is a function of correlation peak height; empirically, we find that the scatter is smallest among 
spectral orders that yield correlation peak heights above a certain threshhold value, typically around 
0.6 (as in the example of Fig. 9b). We thus determine primary radial velocities from each spectrum 
by taking the mean of only those spectral orders with correlation peak heights greater than this 
threshold. Finally, we further refined the primary velocities by iteratively excluding orders that 
deviated by greater than 2a from the mean, which typically resulted in the rejection of at most one 
order per spectrum. 

In Table 5 we report the primary and secondary radial velocities that we derive in this manner. 
In addition to the mean velocities adopted for each star, we give the number of spectral orders 
participating in defining the mean. In addition, we report the O — C velocity residuals compared 
to our best-fit orbit solution (see below) for each of the primary and secondary velocities. The 
uncertainty reported with each velocity is simply the standard deviation of the mean. These 
uncertainties, which average around 0.3 km/s for the primary and 1.4 km/s for the secondary, thus 
represent an estimate of our internal errors only; we discuss instrumental and systematic errors in 
§§4.2.3 and 4.2.4. 

Some of our observations were taken at orbital phases near eclipse, resulting in blended cross- 
correlation peaks that were difficult to separate securely. Indeed, the velocity centroids in these 
cases can be adversely affected for at least a couple of reasons. Cross-correlation peaks in close 
proximity to one another can suffer from "peak pulling" , whereby blending of the two components 
can cause apparent shifts in the peak centroids. In addition, partial blocking of the stellar disc 
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of one star by the other causes the stellar line profiles of the blocked star to become asymmetric, 
narrower, and shifted with respect to the center-of-mass velocity. For our orbit solution analysis we 
thus discarded measurements obtained within 0.05 phase of either eclipse, thereby ensuring that 
the correlation peaks would be reasonably well separated and that occultation effects would not be 
a concern. The only exception to this is the spectrum obtained on UT 2002 Feb 21, which occurs at 
the phase where the secondary is almost completely eclipsed by the primary, so that the secondary 
spectrum is practically not present in the observed light. 

The presence of spots on the stellar surfaces can also introduce distortions that affect the 
observed radial velocities. The importance of this effect for our observations can be evaluated from 
the light curve solutions discussed below, which include the effects of spots. Our model fits to the 
light curves obtained in January 2001 show that the effect is rather small. The theoretical radial- 
velocity distortions due both to stellar occultation (the Schlesinger-Rossiter effect; Rossiter (1924); 
Schlesinger (1909)) and the presence of surface spots are shown in Fig. 10. For the primary star, 
the radial velocity distortions due to spots are undetectable at the precision of our study, having 
an amplitude of at most 0.2 km/s, and less at most orbital phases. For the secondary star, the 
distortions are more pronounced, but still have an amplitude no larger than 1 km/s, and only at 
certain phases. As the light curve data were not obtained simultaneously with the radial velocity 
data, we do not correct the observed velocities for these distortions, noting again that the effect is 
in any case probably negligible. 



4-2.2. Orbit solution 

We determined an orbit solution to the radial- velocity measurements in Table 5 by applying a 
least-squares fit simultaneously to the primary and secondary velocities. The orbit solution involves 
seven parameters: the mass ratio, q; the orbital period, P; the semi-amplitudes of the primary and 
secondary velocities, Ka and Kb] the eccentricity, e; the center-of-mass velocity, 7; and the time 
of primary minimum, To. 

With the system ephemeris determined precisely from the eclipse timing data (see §4.1), we 
chose to fix the values of P and To to those determined photometrically, and forced the orbit to be 
exactly circular (i.e. e = 0). The assumption of circularity is corroborated by the eclipse timings 
analysis (§4.1), and is certainly expected for a short-period binary such as this. In any case, we 
did try a solution with e as a free parameter, but found a best-fit value that was statistically 
insignificantly different from zero. 

Using all of the radial velocities resulted in a fit with r.m.s. residuals for the primary and 
secondary of a a = 0.66 km/s and ctb = 2.55 km/s, respectively. However, the measurement of 
the secondary's velocity from the UT 2001 Dec 18 spectrum showed a residual of —5.75 km/s, 
considerably larger than the other measurements. Indeed, a close look at the distribution of the 
secondary residuals revealed this measurement to be a clear outlier; excluding this measurement 
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yielded a fit with <tb = 1-41 km/s, consistent with the internal errors of the measurements. 

We carefully re-examined (and even reprocessed) the spectrum of UT 2001 Dec 18 and the 
cross-correlation functions to see if we could discern any reasons for the large discrepancy of this 
measurement. The primary's measurement is not obviously discrepant, and neither is the velocity 
of the standard star observed on this night (see §4.2.4). We note however that this spectrum does 
have the poorest S/N of all our HET HRS observations (Table 5). Interestingly, our light curve 
analysis (§4.4) suggests that spots might be significant near the orbital phase of this measurement 
(see Fig. 10). 

In Table 6 we present the orbital parameters of VI 174 Ori resulting from our best-fit orbit 
solution, and display this solution along with the radial- velocity measurements in Fig. 11. 

4-2.3. Instrumental stability 

Our spectroscopic observations allow us to determine an orbit solution with high precision, 
that in turn results in very precise (~ 1%) determination of the stellar masses (Table 6). As 
noted above, the residuals of the radial velocity measurements about the best-fit orbit solution are 
oa = 0.66 km s _1 and <tb = 1-41 km s _1 for the primary and secondary, respectively. Evidently, 
the instrumental stability of the HET HRS is at least as good as 0.65 km s -1 . 

To assess the instrumental stability more carefully, we analyzed our observations of the radial- 
velocity standard star HD 26162, which was observed along with V1174 Ori, in the same way as 
our V1174 Ori observations. Specifically, we cross-correlated each HD 26162 spectrum against the 
same radial-velocity template (HD 18884) that we used to derive the V1174 Ori velocities. The 
results are tabulated in Table 7 and shown graphically in Fig. 12. 

The velocity deviations of HD 26162 relative to its published velocity of 24.8 km s _1 (Udry et 
al. 1999) reveal a high level of stability of the HET HRS system over the 90-day timescale spanned 
by our observations. As Fig. 12 shows, the system does exhibit low-level, secular drifts, but the 
dispersion of these drifts is small. The deviations have an rms of 0.52 km s _1 when we consider all 
of the nights on which VI 174 Ori was observed, and an rms of 0.26 km s _1 when we consider only 
those nights used in our orbit solution 12 . 

This rms of 0.26 km s _1 is only about 1/2 of the rms residuals in our orbit solution fit to the 
primary velocities (0.65 km s _1 ), and about 1/5 of the secondary (1.41 km s _1 ). Evidently, some 
other source of measurement uncertainty dominates the errors in the radial velocity measurements 
of V1174 Ori. These could include: the lower S/N of the V1174 Ori spectra as compared to the 
HD 26162 spectra (Table 3), particularly for the secondary star; the double-peaked nature of the 



The lower rms of this subset of the radial-velocity measurements could be due to chance, given the apparently 
non-random nature of the deviations shown in Fig. 12. 
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Table 5. Radial velocities of V1174 Ori 



# a 


UT Date 


HJD b 


Phase 


R.V.j, 
(km/s) 


(O-C)p 
(km/s) 


R.V. S 
(km/s) 


(0-C) s 

(km/s) 


N P 


N s 


2 


2001 Nov 23 


2236.84827 


0.427 


-9.87 ±0.44 


-0.14 


71.48 ± 1.56 


-2.16 


11 


2 


3 


2001 Dec 10 


2253.80559 


0.863 


83.44 ±0.27 


0.02 






12 





4 


2001 Dec 18 


2261.78887 


0.893 


72.91 ±0.19 


0.04 


-46.57 ± 1.23 


-5.75 c 


11 


3 


6 


2001 Dec 22 


2265.76822 


0.404 


-19.90 ±0.32 


-0.14 


88.17 ± 1.29 


0.68 


15 


6 


7 


2001 Dec 23 


2266.76687 


0.783 


101.81 ±0.34 


0.73 


-80.36 ±1.29 


-0.91 


15 


7 


8 


2002 Jan 26 


2300.67336 


0.652 


88.81 ±0.22 


-0.57 


-62.39 ± 1.43 


0.89 


17 


2 


10 


2002 Jan 28 


2302.66839 


0.409 


-16.27 ±0.23 


1.31 


86.28 ±1.91 


1.81 


20 


3 


11 


2002 Feb 07 


2312.64124 


0.194 


-48.06 ± 0.24 


-0.87 


126.16 ± 1.29 


0.77 


17 


3 


12 


2002 Feb 09 


2314.63298 


0.950 


48.14 ±0.30 


-0.73 


-5.87 ±0.98 


1.45 


17 


2 


14 


2002 Feb 21 


2326.61211 


0.497 


23.23 ±0.29 


-0.34 






20 





15 


2002 Feb 23 


2328.60823 


0.254 


-52.17 ±0.27 


0.30 


131.84 ±1.36 


-0.85 


17 


8 



a See Table 3. 

b Heliocentric Julian Date (2450000±) 

c Measurement discarded in preferred solution (see text). 



Table 6. Results of VI 174 Ori orbit solution 



Parameter 


Value 


P(d) 


2.634727 ± 0.000004 


a a (km s r ) 


0.66 


og (km s _1 ) 


1.41 


1 


0.7238 ± 0.0057 


e 




7 (km s _1 ) 


25.27 ± 0.19 


K A (km s- 1 ) 


77.75 ± 0.29 


K B (km s- 1 ) 


107.41 ± 0.67 


asini (R©) 


9.638 ± 0.041 


M A sin 3 i (Mq) 


1.005 ± 0.015 


M s sin 3 i (Mq) 


0.728 ± 0.008 



a Adopted. 
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V1174 Ori cross-correlation functions; differences in the track length of the V1174 Ori observations 
as compared to the HD 26162 observations (1800 sec vs. 10 sec); a source of astrophysical noise (e.g. 
surface activity such as spots, flares, etc.); and/or the presence of an unseen tertiary companion (see 
§5.1). In any case, instrumental stability does not appear to be the dominant source of uncertainty 
in our radial velocity measurements. 

4-2.4- Systematic errors 

For this experiment, where the primary goal is to determine accurate stellar dimensions for 
comparison to stellar evolutionary models, it is critical that we ascertain the extent to which 
systematic effects may limit the accuracy of the absolute stellar dimensions that we determine. 

The semi-amplitudes of the primary and secondary velocities [Ka and Kb) are the parameters 
from the spectroscopic orbit solution that directly determine the stellar masses, so the masses are 
directly sensitive to systematics in these quantities. One possible source of systematic error in these 
quantities is spectral-type mismatch between V1174 Ori and the radial-velocity template. 

There is in fact some evidence for a systematic offset of the secondary velocities relative to 
the primary. As Fig. 12a shows, the velocities of the radial-velocity standard HD 26162 show a 
median velocity in our observations that is 0.36 km/s larger than the published value. We speculate 
that this offset is due to the fact that HD 26162 has an earlier spectral type (SpT = K2 III) than 
the radial- velocity template star (SpT = M2 III). Our best estimate for the spectral types of 
the components of V1174 Ori are K4.5 and Ml. 5 (see §4.3.3), so a similar systematic shift might 
be expected between the measured velocities for the two components of V1174 Ori, leading to 
systematic errors in the radial velocity semi-amplitudes. 

We tried an orbit solution where we allowed for an arbitrary systematic offset between the 
primary and secondary velocities to be fit. We find a best-fit value of 0.35 =L 0.55 km s _1 for the 
offset, which is interesting in its similarity to the offset of 0.36 km s _1 we find above. Including 
this free parameter and its uncertainty results in slightly different stellar masses than those listed 
in Table 6, and with slightly larger uncertainties: Ma sin 3 i = 1.004 ± 0.016 M and Mb sin 3 i = 
0.727 ± 0.009 M Q . Thus, while our orbit solution formally allows for an offset of order 1 km s _1 at 
the 2(7 level, it is not statistically significant and its effect on the final uncertainties on the derived 
stellar masses is very small. For our subsequent analyses and discussion, we adopt the orbit solution 
provided in Table 6. 

4.3. Spectral analysis 

Having determined precise dynamical masses from our spectroscopic orbit solution, we now 
require stellar radii and effective temperatures for a complete analysis of the VI 174 Ori system 
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and for comparison with PMS stellar evolutionary tracks. Our light curve analysis (§4.4) provides 
precise values for the radii, the inclination, and the ratio of effective temperatures. Thus we need 
to establish the effective temperature of one of the components independently. 

In this section we conduct a more detailed spectral analysis to determine the primary star's 
effective temperature. To begin, we measure the rotational velocities of the two stars and determine 
their flux ratio. This information guides our spectral classification of the primary by establishing 
the rotational broadening of the primary spectrum and the extent to which the secondary spectrum 
contaminates the lines that we use to classify the primary. The rotational velocity and flux ratio 
information will also prove useful in constraining our light curve analyses below. Finally, we measure 
the strength of two spectral features of interest in studies of young stars, Li and Ha, the former of 
which serves to establish stellar youth and the latter of which serves as a proxy for chromospheric 
activity and/or accretion. 

4-3.1. Rotational velocities 

We measure the rotational velocities, v sini, of the two components of V1174 Ori in two ways. 
First, the cross-correlation functions we used to measure radial velocities (§4.2.1) contain rotational 
information in the sense that the peak widths are related to vsini. Second, we use the ratio of the 
widths of the primary and secondary Li lines to constrain the ratio of v sin i of the two components. 

Because the widths of the cross-correlation peaks reflect the convolution of the intrinsic stellar 
vsini with both the instrumental response function and the radial- velocity template spectrum, 
the peak widths do not directly yield a measure of vsini. We calibrated the width of the cross- 
correlation peak to v sin i by artificially rotationally broadening our narrow-lined (lines unresolved) 
spectral standards (see Table 3) and then cross-correlating these broadened spectra against the 
same radial-velocity template used with our V1174 Ori spectra. This was done separately for each 
spectral order. We find that our cross-correlation functions permit us to discern vsini values of 
> 10 km/s. 

Using this calibration, we measured vsini values for the components of V1174 Ori from each 
spectrum, averaging together the values from the multiple orders in a given spectrum. The resulting 
vsini values are 26.2 ± 0.6 km/s and 20.5 =L 1.6 km/s for the primary and secondary, respectively. 
Here the uncertainties are the standard deviation of the mean. 

We wish to emphasize that the uncertainties quoted above are strictly internal. Furthermore, 
we caution that we have not used vsini standards for an external calibration. For our present 
purposes, what is robust here is that the rotational velocities for the stars are ~ 20-25 km/s, and 
that the primary star rotates more rapidly than the secondary. This latter point is corroborated 
by our analysis of the primary and secondary Li lines below (§4.3.4). 
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4-3.2. Flux ratio 

We use the temperature-insensitive Ca II lines at AA6102,6122 to measure the flux ratio between 
the primary and secondary components of V1174 Ori. We fit the V1174 Ori spectrum of UT 2002 
Feb 07 by adding together at various ratios the spectra of our K3 and K7 standards 13 (see Table 3) 
at the appropriate radial velocities and with the appropriate rotational broadening. We chose the 
UT 2002 Feb 07 spectrum because on this date the secondary was shifted a full 174 km/s (3.5A at 
this spectral region) redward of the primary. The two components are thus well separated, and the 
A6102 line of the secondary is shifted into a region of the spectrum (6105. 5A) that is free of lines 
from the primary. 

We find a best fit to the Ca II lines in the V1174 Ori spectrum when we combine the two 
templates in a 6:1 ratio, (Fig. 14), although ratios of 5:1 and 7:1 also fit the observed spectrum 
reasonably well. 

4-3.3. Spectral classification 

Our high-resolution HET HRS spectra permit spectral classification of the primary from careful 
examination of temperature-sensitive line ratios. As demonstrated by Basri & Batalha (1990) and 
by Lee (1992), the ratio of Ni I A6108 to V I A6112, the ratio of V I A6040 to Fe I A6042, the ratio of 
V I A6058 to Fe I A6056, and the ratio of Sc I A6210 to Fe I A6200 provide a good determination of 
spectral class from high-resolution spectra. These line ratios were also used by Steffen et al. (2001) 
in their analysis of the PMS binary NTT 045251+3016. 

We have carefully examined these line ratios in our spectra that have the secondary sufficiently 
Doppler-shifted with respect to the primary. We find a spectral type for the primary of K4-K5, 
with an uncertainty of ~l/2 subclass. As an example, in Fig. 13 we compare the first line pair 
above from our HET HRS spectrum of UT 2002 Feb 07 with four spectral standards of type Kl, 
K3, K4, and K5. The observed line ratio indicates a K4.5 spectral type, but after correcting the 
primary's V I line for contamination by the secondary's Ni I line (using the flux ratio determined 
above), we find a K4 spectral type in this spectrum. 

Taking all of these line-pair ratios from all of our observations into account, we find a mean 
spectral type for the primary of K4.5±0.1. Due to the quantized nature of the spectral classification 
process (0.5 sub-type), this uncertainty is likely unreastically small, so we adopt a more conservative 
uncertainty of 0.5 spectral sub-type. Using the SpT-T e g- conversion of Schmidt-Kaler (1982), this 
implies an effective temperature for the primary of T e fj = 4470 ± 120 K. 



Our use of the K3 and K7 standards is based on the preliminary spectral types that we determined for the two 
components of VI 174 Ori, but as the Ca lines used are temperature-insensitive any two standards of roughly mid-K 
to early-M spectral type would suffice. 



-19- 



Using this together with the ratio of T e s determined from our light curve analysis below (§4.4), 
we infer that the secondary has an effective temperature of T e g- = 3615 K, implying a spectral type 
of ~M1.5. 

4-3.4- Lithium 

The A6708 line of Li is a commonly used indicator of stellar youth, as stars rapidly destroy 
their natal Li content during the first ~fewxl0 7 yr of their evolution. For very young stars, the 
Li line at 6708A can be very strong, with equivalent widths from ~ 300 mA for early K stars to 
~ 750 mA for mid-M stars (Dolan & Mathieu 1999). Indeed, in our HET HRS spectra, the Li line 
is one of the few lines from the secondary that can be unambiguously identified by eye. 

We have measured the equivalent widths of the primary and secondary Li lines in those spectra 
where the two stars are sufficiently Doppler-shifted relative to one another. As an example, the 
spectrum of UT 2002 Jan 28 has the secondary at a radial velocity of +102.55 km/s relative to 
the primary. The Li equivalent width (EQW) we measure for the primary star in this spectrum is 
EQW(Li) A = 398 mA. We measure the Li EQW of the secondary star to be EQW(Li) B = 78 mA 
after correcting for the contribution of spectral features from the primary. These EQW measures 
also include a 5% correction for the presence of third light, based on our light curve analysis (§4.4.3) 
which indicates 0% third light at V and 10% at /. 

If intrinsically equal, these EQW values would imply a flux ratio of 5.1. However, adopting 
spectral types of K4.5 and Ml. 5 for the two components implies that the intrinsic EQW of the 
secondary's Li line is ~ 25% greater than that of the primary (Dolan & Mathieu 1999). Correcting 
the secondary's EQW(Li) accordingly results in a flux ratio of 6.3, consistent with that derived 
from our analysis of the temperature-insensitive lines above (§4.3.2). Thus from this spectrum we 
infer intrinsic equivalent widths for the two stars of EQW(L1)a ~ 465 mA and EQW(Li)e ~ 545 
mA. The other spectra for which we could cleanly separate the primary and secondary Li features 
yield values consistent with these. Adopting a flux ratio of 6:1, we find from the ensemble of our 
spectra EQW(Li) A = 480 ± 14 mA and EQW(Li) B = 537 ± 26 mA. 

We can use these EQW measurements to infer Li abundances for the two stars. From the 
Li curve-of-growth analysis of Pavlenko & Magazzu (1996), we find N(Li)a = 3.08 ± 0.04 and 
N(Li)b = 2.20 ± 0.06 assuming LTE, and N(Li) A = 3.00 ± 0.05 and N(U) B = 1.98 ± 0.09 for the 
non-LTE case. Note that these abundance estimates are likely to be slightly lower than the actual 
abundances, as the curve-of-growth values used assume logg = 4.5, whereas the components of 
V1174 Ori have logg « 4.2 (see §5.1). In addition, the uncertainties listed reflect only the error in 
the measurement of the EQW; including the uncertainty in the stellar T e g increases the abundance 
uncertainties to 0.2 dex. 

Thus, the primary star's Li abundance is consistent with at most 0.2-0.3 dex of depletion 
(assuming a "cosmic" Li abundance of logiV(Li)= 3.3), whereas the secondary shows > 1 dex of 
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depletion. As discussed by D'Antona & Montalban (2003) the pattern of PMS Li depletion as a 
function of stellar mass and age can provide important constraints on PMS tracks. We will return 
to the issue of Li abundances in §6.5. 

Finally, we can also use the primary and secondary Li lines to infer the relative rotation rates 
of the two stars. We have already derived v sin i values from the widths of the cross-correlation 
peaks (§4.3.1). Using the ratio of the widths of the Li lines, we can independently infer the ratio of 
vs'mi. We have measured this ratio from our spectra that have the primary and secondary Li lines 
sufficiently well separated to allow a clean measurement of their FWHM, and find a vsini ratio of 
1.21 ± 0.04, consistent with the vsini values determined from the cross-correlation analysis. 

4.3.5. Ha 

Emission of the Ha line is also of interest in studies of young stars. Indeed, emission stronger 
than EQW(Ha) m 10 A is a defining criterion for classical TTS, which are believed to be PMS 
stars actively accreting from circumstellar disks. Even non-accreting TTS often show some amount 
of Ha in emission, which is thought to indicate the presence of strong chromospheric activity. 

Based on the strength of Ha emission, VI 174 Ori does not appear to be an actively accreting 
system, with neither component exhibiting Ha emission in excess of a few A EQW in any of our 
spectra. In Fig. 15, we show the spectrum of UT 2002 Feb 23 in the vicinity of Ha. This spectrum 
has the components of VI 174 Ori at their greatest radial-velocity separation, making it possible 
to detect their individual Ha emission apart from the very strong nebular emission. The observed 
equivalent widths measure EQW(Ha)A = — 0.4A and EQW(Ha)e = —4.2 A, after correcting for 
the continuum of the companion. 

4.4. Light curve synthesis and modeling 

Having determined an ephemeris for V1174 Ori, and having conducted a thorough spectral 
analysis of the system including determination of a double-lined orbit solution, we now analyze 
our multi-epoch, multi-band light curves (Table 1, Figs. 2-7) in detail to determine the orbital 
inclination of the system, T c g- for the secondary, and absolute radii for both components. Our 
analysis procedure makes use of a Wilson-Devinney (WD) based algorithm that models the observed 
photometric variations with synthetic light curves that include the effects of limb darkening, gravity 
brightening, mutual illumination and occultation, and spots. In addition, the code permits the use 
of either simple Planckian fluxes or detailed model atmospheres. Along the way, we use information 
provided by our spectroscopic analyses (vsini, flux ratio, primary T e g) to further constrain the 
model light curve parameters. 
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4- 4-1- Initial considerations 

Although the light curves shown in Figs. 2, 4 and 6 are typical for a system with non-distorted 
components, apparently well-detached from their Roche lobes, they also reveal the presence of 
variations not related to the eclipses. From careful inspection of the light curves it becomes evident 
that these variations are most likely due to star spots that change from one epoch to the next, 
hindering the use of our photometric data from all epochs in a single, integrated analysis. 

Thus, to begin, we used the data from the epoch with the most complete set of light curves — 
those obtained in Jan 2001 in 5 filters (vbyVI) — to produce an initial set of light curve solutions. 
These initial light curve solutions were calculated without spots by fitting simultaneously the 5 
light curves from Jan 2001 together with the radial velocity data from Table 5. 

Our analysis was performed with the WD light-curve synthesis program (Wilson & Devinney 
1971; Wilson 1993a,b). Besides other modifications (Vaz et al. 1995; Casey et al. 1997), our version 
of the WD code was implemented with the possibility of using the SIMPLEX solution method 
(Kallrath & Linnel 1987; Vieira 2003) as an alternative to the traditional least-squares method. 
Limb-darkening coefficients for both eclipsing components were described using the data from 
van Hamme (1993). All three limb-darkening laws (linear, linear-log, and linear-square root) were 
examined, with the coefficients calculated by bi-linear interpolations using the current values of log g 
and T e ff. The linear law gave marginally better solutions and was therefore adopted. Bolometric 
albedoes for both components were set equal to 0.5 as appropriate for atmospheres in convective 
equilibrium. Gravity-brightening exponents were calculated using the local value of T c g for each 
point on the stellar surfaces, taking into account mutual illumination, according to Alencar & Vaz 
(1997) and Alencar et al. (1999). The radiated flux of both stars is described using the atmosphere 
model tables of Buser & Kurucz (1992), updated relative to the original versions of the WD code. 
As the T e ff of the secondary is lower than the minimum T e g included in the published tables (4000 
K), we extrapolated linearly in a log(theoretical magnitude) vs. log(temperature) scale. Both the 
radial velocities and the photometric times of minima indicated no sign of orbital eccentricity, so 
we assumed circular orbits in the light curve analysis. 

In performing these initial calculations, we applied least-squares differential corrections be- 
tween successive iterations to the orbital inclination, the secondary T e fj, the primary luminosity, 
both stellar surface gravitational pseudopotentials (stellar radii), the center-of-mass radial velocity, 
the mass ratio and an arbitrary phase shift to be added to the observations. The luminosity of 
the secondary was calculated from its size and T c fj. In addition, the limb-darkening coefficients, 
normalization magnitudes, surface gravities, and individual velocity amplitudes were all updated 
between successive runs to correspond to the solution from the previous iteration. 

Initially, the WD code settled on a solution in which the secondary star was larger than the 
primary. However, this solution violated the observed luminosity ratio (§4.3.2) as well as the 
observed rotational velocities if synchronous rotation is assumed. By forcing the luminosity ratio 
constraint, the code found another solution in which the primary is larger. This solution matches 
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the observed luminosity ratio and predicts rotational velocities similar to those observed. This 
solution thus served as a starting point for introducing additional parameters into the analysis, 
such as spots, as we now describe. 

4- 4- 2- Star spots 

After determining an initial solution as described above, we introduced spots in the WD 
calculation in order to model the out-of-eclipse variations observed. We started by introducing one 
cold spot on each component. For each spot on each star, four parameters are added to the WD 
calculation: the longitude and co-latitude of the spot, the spot radius (as seen from the center of 
the star), and the spot temperature factor (< 1 for cold spots and > 1 for hot spots). 

In this phase of the calculation, we found the SIMPLEX method to be much more efficient 
in finding best-fit solutions than the least-squares method, the latter being very sensitive to small 
corrections to the spot parameters, leading to frequent divergence of the solution. Although the 
WD algorithm can deal with spots that move in longitude across the surface of the star, we assumed 
that the spots rotate with the same angular velocity as the star. After convergence of each set of 
spots, we tried many different combinations of spot number and configurations, trying to minimize 
the number of spots on each star. Note that the WD code treats regions in which multiple spots 
overlap by multiplying their temperature factors 14 . 

Ultimately, both "cold" and "hot" spots were needed on each star to reproduce the light 
curves of Fig. 2. The cold spots made it necessary to extrapolate in the atmosphere tables, as 
mentioned above for the secondary star, for some points on the surface of the primary star as well. 
The calculations with atmosphere tables were compared with calculations using blackbodies, and 
we concluded that the extrapolation done (linearly in a logarithmic scale) was reasonable. Thus, 
except for the JHK light curves for which Buser & Kurucz (1992) did not publish atmosphere 
calculations, we thus chose to use atmosphere models instead of simple Planck functions. 

There exist multiple combinations of spots that can reproduce the observed out-of-eclipse 
variations reasonably well. We believe that the 3 spots on each star used in our final solution 
(Table 9) are physically realistic. Moreover, our use of light curves in five different bands makes 
the spot solution more robust, as we insisted that one spot configuration reproduce all of the light 
curves simultaneously. In any case, the goal of including spots in the light curve modeling is not 
to determine the properties of the spots per se, but rather to better constrain the morphology of 
the light curves around the eclipses and thereby to minimize the uncertainty in the parameters of 
true import: the system inclination, the ratio of T e g, and the stellar radii. 



14 This differs from, e.g., the Nightfall program, which uses the average of the temperature factors. See 
http://www.lsw.uni-heidelberg.de/users/rwichman/Nightfall.html . 
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44-3. Third light 

The Jan 2001 light curves proved impossible to fit together unless we included ~10% third 
light in /. The V light curve required a slightly negative amount of third light, which is physically 
meaningless. We thus forced L%{V) = 0, but left this parameter free to be adjusted at the other 
wavelengths. The final solution is shown in Figs. 2 and 3 and in Tables 8 and 9. With the exception 
of /, the significance of the L3 values found for the light curves is questionable as they do not follow 
a Planckian distribution with wavelength. They are probably artifacts of the data-reduction process 
(e.g. background subtraction in the aperture photometry, scattered light, etc.). 

4-4 -4- Final solutions 

In order to test the validity of the stellar parameters found above from the light curve solutions 
for the Jan 2001 data, we applied the same analysis procedure to the other sets of light curves (Table 
1), now forcing all the parameters, except those corresponding to the star spots, to remain fixed. 
The resulting fits for the / light curves obtained in Nov 2001 and the BVI light curves obtained in 
Nov-Dec 2002 are shown in Figs. 4 and 5. The characteristics of the spots used are given in Table 
9. The amount of third light at the different wavelengths were left as free parameters. L^(V) and 
Lz(B) tended to become slightly negative, while L^{I) was once again ~ 10%. We thus used no 
third light for V and B, and adopted the value for L^(I) shown in Table 8. 

The /-band discovery light curve obtained in Dec 1994 (Table 1) was also fit with the solution 
of Table 8 and a different spot configuration. The amount of third light was again fixed at the 
value shown in Table 8. Similarly, the JHK light curves from Carpenter, Hillenbrand, & Skrutskie 
(2001), while not providing complete phase coverage, were also fit with this solution. No third light 
was included in the JHK light curves because the third light parameter mainly affects the depth 
of the minima, and unfortunately the minima were not covered by the Carpenter, Hillenbrand, & 
Skrutskie (2001) observations. The solutions are shown in Figs. 6 and 7, and the spots used are 
listed in Table 9. 

Overall, the quality of the synthetic light curve fits to the data is excellent. Indeed, all 5 bands 
from the Jan 2001 observations are fit extremely well, with only minor systematic effects visible in 
the residuals (Fig. 3). Small systematic residuals do persist in some of the light curves at other 
epochs. For example, in the Nov-Dec 2002 /-band data there are 2 points that deviate by more 
than 0.02 mag relative to the theoretical light curve near secondary minimum (Figs. 4 and 5). The 
data points close to secondary minimum in V and B also show small deviations, although at only 
~ 0.01 mag these are not by themselves statistically significant. A few deviant data points also 
appear at other phases (e.g. phase 0.3 in I). These might be indicative of erratic variations as have 
been reported in, e.g., AK Sco (Alencar et al. 2003) and TY CrA (Casey et al. 1997). Indeed, the 
most discrepant point near phase 0.5 in the Nov-Dec 2002 I-band light curve was obtained by itself 
more than two weeks later than the other points near that phase; perhaps there was small amount 
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of spot evolution or some other small change in the system over the course of the Nov-Dec 2002 
epoch. 

A limitation in our ability to analyze small-scale photometric variability over long timescales 
is that our differential photometric observations were not all reduced in one integrated procedure, 
but separately for each epoch. The Honeycutt (1992) algorithm we employ determines the dif- 
ferential photometric variations relative to a statistically defined ensemble of comparison stars, 
and this statistical ensemble is not necessarily common to all epochs. Consequently, the photo- 
metric calibrations from one epoch to the next are not absolutely tied to one another. In both 
AK Sco and TY CrA, secular changes in system brightness from one epoch to the next have been 
attributed to variable obscuration by dust. Unfortunately, we cannot determine from our current 
analysis whether such global changes in system brightness are taking place in VI 174 Ori over long 
timescales. 

Even so, the light curve analyses presented here show clear evidence for changes with time in 
both the number and the characteristics of the spots on both components. The time coverage of 
the observations was not adequate to follow the evolution of the spots in detail, but there are hints 
in Table 9 that some of the spots may persist on month timescales, drifting across the surfaces of 
the stars. The hot spots apparently come and go more sporadically and, besides occurring more 
frequently on the secondary than on the primary, they tend to be smaller than the cold spots. 
Finally, although the rotation rates of Table 8 indicate orbital synchronism within the errors, one 
of the reasons for spot drifting may be that the stars are not yet completely synchonized. The stars 
may also present differential rotation, as in the Sun and as recently determined for the rapidly- 
rotating RS CVn binary UZ Librae (Olah et al. 2003). 

In the end, while some questions remain as to the origin of the variability observed in some of 
our light curves at the finest level of detail, the overall match between the synthetic light curves 
and the data is excellent (Figs. 2-7), and for the remainder of our analysis we adopt the solution 
listed in Table 8. 

5. Results 

5.1. Fundamental stellar parameters 

Combining our spectroscopic and photometric analyses presented above, we can establish all of 
the fundamental stellar properties for both components of VI 174 Ori. We report these parameters 
in Table 10. 

We can in principle determine the extinction, Ay, toward the system by fitting the observed 
colors (Table 2), assuming main-sequence intrinsic colors (Kenyon & Hartmann 1995) and bolo- 
metric corrections (Popper 1980) for the primary and secondary components of V1174 Ori. There 
is not a single value of Ay that simultaneously reproduces all of the observed colors. Formally, 
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however, we find a best-fit value of Ay = 0.32, which results in a distance to the system of 419 pc 
when we normalize to the observed V magnitude (see §5.2 below). The colors calculated with this 
extinction and distance are compared to the observed values in Table 2. 

The agreement between the observed and calculated colors is best in B — V, where the observed 
flux is dominated by that of the primary star, and in J — H and H — K, where the colors become 
largely degenerate for dwarfs of mid-K to early-M spectral type. The observed U — V color is bluer 
than predicted by 0.11 mag, perhaps the result of a mild UV excess from one or both of the stars, as 
is common among chromospherically active PMS stars, and possibly consistent with the observed 
Ha emission (see §4.3.5) and hot spots (§4.4.2). 

The observed V — I color is somewhat redder than predicted (by 0.14 mag), roughly consistent 
with the level of third light inferred from our light curve modeling (§4.4) in the I band (~ 10%). 
Similarly, the observed J magnitude is 0.30 mag brighter than predicted. 

We can improve slightly the fit to the observed V — I color and J magnitude by increasing the 
Ay. For example, adopting Ay = 0.42 would predict V — I = 1.52 and J = 11.46. Of course, this 
worsens the agreement in the other colors, most notably in U — V, and results in a distance (401 
pc) less compatible with that typically assumed for the ONC (~ 470 pc). 

Together, the V — I and J discrepancies suggest the presence of a third, low-mass star in the 
system that contributes flux noticeably only at A > 0.9/im. This would imply that V1174 Ori is in 
fact a hierarchical triple system, with the tertiary in a sufficiently wide orbit so as to not manifest 
itself in the observed radial velocities (§4.2.1). Curiously, TY CrA (Casey et al. 1998), one of the 
other three PMS eclipsing binaries in the literature, is also reported to be a triple system. 

5.2. Distance to V1174 Ori 

Since we measure the luminosities of the components of V1174 Ori directly from T c g and R, 
we can infer a distance to the system assuming a particular Ay derived from the system colors 
(see above). The distance we derive to VI 174 Ori of 419 ± 21 pc differs somewhat from that which 
has been quoted for the ONC in the recent literature. For example, Hillenbrand (1997) adopt a 
distance of 470 pc, and Bally, O'Dell, & McCaughrean (2000) use a distance of 460 pc. 

Various other authors have determined distances to the ONC and to the various subgroupings 
of Orion OBI, ranging from 400-500 pc. Based on photometry of several tens of bright O and B 
stars, Warren & Hesser (1978) find a distance of 480 pc to the Trapezium cluster (Ori OBld) and 
430 pc to the foreground Ori OBlc, while Anthony-Twarog (1982), using the same data, find a 
distance of 434 pc for Ori OBld. Using different photometric indices, Brown, de Geus, & de Zeeuw 
(1994) find a distance of only 400 pc for Ori OBlc. Perhaps the most robust constraint is that 
provided by the BN-KL masers, situated just behind the ONC, at 480 pc (Genzel et al. 1981). 

Given precisely determined luminosities for the two eclipsing components of V1174 Ori, our 
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distance determination depends primarly upon Ay, with lower values placing the system farther 
from the Sun. Based on our analysis of the observed colors, we do not believe that Ay can be 
much lower than what we report above. Using the value of Ay reported in Table 2, the discrepancy 
between the observed and predicted V — I color is already consistent with the amount of third 
light that we find for the I band from our light curve analysis. Decreasing Ay would increase the 
discrepancy in V — I, which would be incompatible with our light curve solution. 

Once the value of Ay is determined, we compute the distance via the distance modulus in V. 
Our U-band measurement has a formal uncertainty of 0.03 mag, which translates into an additional 
distance uncertainty of 6 pc. Of course, this uncertainty does not include any possible influence of 
spots which, from our analysis above (§4.4.2) has been observed to be ~ 0.1 mag peak-to-peak in 
V. Our calibrated photometry was obtained in 1995, an epoch at which we unfortunately do not 
have light curves that could inform us as to the influence of spots on the calibrated magnitudes. 
If, for example, a cool spot were acting to depress the F-band flux by 0.1 mag, then our derived 
distance would be underestimated by 20 pc. 

Another, if subtler, consideration in our distance determination is the value assumed for the 
solar bolometric magnitude, Mb i,©- Lower values of Mboi,© (i.e. brighter) result in VI 174 Ori 
being farther from the Sun. We have adopted the value of Mbd,© = 4.59 used by Caloi, D'Antona, 
& Mazzitelli (1997), which produces a distance roughly consistent with the other distance determi- 
nations discussed above. However, the value currently adopted by the IAU, Mboi,© = 4.75, results 
in a distance to V1174 Ori of 390 pc. 

Based on this discussion, we consider it likely that V1174 Ori is a member of the older, 
foreground Ori OBlc at a distance of ~ 400 pc, consistent with the slightly older age (~ 10 Myr) 
that we find below. Indeed, the position of V1174 Ori approximately 0.3 degrees to the south of 
Q 1 C Ori places it just outside the region defined as Ori OBld by Blaauw (1964) and in the region 
defined as Ori OBlc. 

6. Comparison to theoretical pre main-sequence evolutionary tracks 

In this section, we use the empirically determined masses, radii, effective temperatures, and 
luminosities of the primary and secondary components of V1174 Ori to test the predictions of PMS 
stellar evolutionary models. We consider in turn four sets of PMS tracks commonly used in the 
literature: Baraffe et al. (1998) (§6.1), Palla & Stahler (1999) (§6.2), Siess, Dufour, & Forestini 
(2000) (§6.3), and Montalban et al. (2003) (more recent version of the D'Antona & Mazzitelli (1997) 
tracks; §6.4). 

These models differ primarily in their choice of atmospheres/opacities, the equation of state, 
and in their treatment of convection. With respect to the former, the models essentially differ in 
terms of whether the atmosphere models used are gray or non-gray, and whether molecular opacity 
effects are included at low effective temperatures. The differences with respect to convection are 
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somewhat more subtle. Except for one set of models that uses the Full Spectrum Turbulence 
(FST) formalism of Canuto et al. (1996), the models we consider in this paper use a mixing-length 
theory (MLT) approach, typically characterized by the mixing parameter, a, defined as the ratio 
of the convective mixing length to the pressure scale-height, i.e. a = l m i x /H p . But as discussed by 
Montalban et al. (2003) and D'Antona & Montalban (2003), to fully characterize an MLT model 
requires specification of a larger set of ingredients than just the single parameter a. These are: 
the mixing parameter used in the atmosphere calculation (a a t m ); the mixing parameter used in 
the calculation of the interior (a- m , which can be set equal to cc a tm); and the optical depth in the 
atmosphere at which these two calculations are matched up (r p h). Requiring a fit to the present- 
day Sun generally requires more efficient convection in the interior than in the atmosphere, that is, 
am > Oatm (Baraffe et al. 1998; Montalban et al. 2003). For a given choice of a atm and a m , larger 
T p h produce lower temperatures at the base of the convection zone, resulting in less Li depletion 
prior to the main sequence. The parameters used by the various models considered below are 
summarized in Table 11. All of the tracks considered assume solar metallicity. 

For each set of tracks, we first consider the positions of the components of V1174 Ori in the 
H-R diagram (T e g vs. log L) and compare the stellar masses inferred from the models to those we 
have measured. We also consider whether the theoretical isochrones lie parallel to the observed 
positions. 

Comparing the observations to the models in the plane of the H-R diagram is not necessarily 
the most robust test of the models, however, and certainly does not take full advantage of the 
accuracy of the measurements. Placing the stars in the H-R diagram requires deriving L from 
R and T e fj. While the measurement errors in the radii are small, the T c g values have not been 
determined directly but rather converted from the observed spectral type (for the primary) and 
then from the T e g ratio determined from the light curve analysis (for the secondary) . Consequently, 
the uncertainties in the derived luminosities are relatively large. Thus for each set of tracks we also 
do the comparison in the mass-radius plane. In this way, not only do we preserve the accuracy of 
the empirical masses and radii, we also test what is perhaps the most fundamental prediction of 
the models: the mass-radius relationship. 

One final comment is in order, regarding the placement of the primary and secondary stars in 
the H-R diagrams shown in Figs. 16-21. In each H-R diagram we plot two nested boxes around 
the measured position of the primary star. The inner box represents the region of the (logL, T e g) 
space allowed by the combination of the uncertainty in the primary star's T e g (0.5 spectral sub- 
type; §4.3.3) and la uncertainty in its radius (0.015 R0; Table 10). The outer box corresponds 
to twice the uncertainty in T e g and 2a uncertainty on the radius. This outer box should thus be 
regarded as the region of high confidence for the location of the primary star in the H-R diagram. 
We plot similar nested boxes about the position of the secondary star, but the range of T e R here 
is determined by keeping the secondary-to-primary T e g ratio fixed at the value determined from 
the light curve analysis (§4.4), the uncertainty of which is minuscule compared to the uncertainty 
in the T e g of the primary star. Because the two stars' positions are linked via the T e g ratio, it 
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is important to bear in mind that the positions of the two stars in the H-R diagram are very 
strongly correlated. While it may be tempting to adjust the position of each star within its own 
bounding box independently, the positions of the two stars relative to one another are in fact rigidly 
constrained. 



6.1. BCAH98 

We begin by considering the PMS tracks of Baraffe et al. (1998). Of the models that we 
consider here, these are perhaps the most generally favored in the recent literature (White et al. 
1999; Simon, Dutrey, & Guilloteau 2000; Luhman et al. 2003). The tracks use the non-gray NextGen 
atmospheres of Hauschildt, Allard, &: Baron (1999) and treat convection using MLT. These models 
are shown in Figs. 16 and 17 for two different values of aj n (1.0 and 1.9, respectively). Note that 
both models use NextGen atmospheres computed with a a t m = 1-0 (see Table 11); thus, these 
models differ in their convection properties only interior to an optical depth of r p h = 100. The 
tracks with a; n = 1.9 are designed to match the properties of the present-day Sun. 

Neither of these Baraffe et al. (1998) models succeed in simultaneously reproducing the ob- 
served positions of both components of V1174 Ori in the H-R diagram. The tracks with = 1.0, 
favored by Steffen et al. (2001) in their analysis of the PMS binary NTTS 045251+3016, predict a 
mass for the V1174 Ori primary of > 1.1 M . The position of the primary excludes the 1.01 M 
track with > 3a confidence. Meanwhile, the 0.73 M track is consistent with V1174 Ori secondary 
at the ~ 1.5a level. 

The tracks with a- ul = 1.9, favored by the recent analysis of Luhman et al. (2003), are possibly 
consistent with the position of the 1.01 M primary of V1174 Ori, but only by pushing the primary 
near the lower extreme of its 2<7 bounding box. Forcing agreement with the primary in this way 
necessarily forces the secondary to the lower extreme of its 2a bounding box as well. This would 
imply a mass for the secondary star of ~ 0.4 M , completely inconsistent with its measured mass 
of 0.73±0.01 M . The 0.73 M track is ~ 400 K too hot relative to the position that the secondary 
would take if the primary were forced into agreement with its mass track. 

Of the BCAH98 models considered here, only the a- m = 1.9 tracks possess an isochrone in the 
H-R diagram (at ~ 7 Myr) consistent with coevality for the components of V1174 Ori. The other 
tracks make the secondary appear over-luminous relative to the primary. 

In the mass-radius (M-R) plane, both sets of tracks imply an age for the V1174 Ori system 
of ~ 10 Myr, though none of the tracks possess an isochrone precisely parallel to the empirical 
isochrone defined by the components of V1174 Ori. Interestingly, this is true of the a = 1.9 tracks 
as well, despite the fact that these tracks yield coevality for the two stars in the H-R diagram. 
Apparently, the empirical measurements and the theoretical isochrones translate differently between 
the M-R plane and the H-R diagram, almost certainly the result of differences in the T e g- scale. 
This is most readily seen in the a- m = 1.0 tracks where, in the M-R plane, the secondary is situated 



-29- 



below the 10 Myr isochrone and the primary is above it; in the H-R diagram their positions relative 
to this isochrone are reversed. 

6.2. PS99 

The PMS tracks of Palla & Stahler (1999) have been found to show reasonably good agreement 
with PMS mass measurements in a number of previous analyses (Steffen et al. 2001; Simon, Dutrey, 
& Guilloteau 2000; Palla & Stahler 2001). These models arguably use the most physically realistic 
initial conditions of any of the models, that of the stellar birthline (Stahler 1983). They use a 
gray atmosphere approximation, and use the Alexander & Ferguson (1994) and Iglesias & Rogers 
(1996) opacities for T e g below and above 10 4 K, respectively. Convection is treated with MLT and 
ctin = 1.5. These tracks are compared with the observations of V1174 Ori in Fig. 18. 

The PS99 models show good agreement with the position of the primary star in the H-R 
diagram, with the 1.01 M track passing within the inner bounding box of the primary's position. 
However, the position of the secondary star excludes the 0.73 M Q track with > 3er confidence. As 
with the BCAH98 tracks, the 0.73 M Q track is ~ 400 K too hot compared to the position of the 
VI 174 Ori secondary. 

The PS99 models moreover yield ages for the two components of V1174 Ori that are non- 
coeval. The tracks give ages of 10 Myr and 5 Myr for the primary and secondary, respectively. 
Interestingly, this non-coevality in the H-R diagram is not borne out in the M-R plane, where the 
components lie precisely parallel to the theoretical isochrones (9 Myr in this case). Indeed, of the 
models considered in this paper, these models show the best agreement with the slope of the M-R 
relationship as defined by the components of V1174 Ori. 

6.3. SDFOO 

Next, we consider the PMS models of Siess, Dufour, & Forestini (2000). These models, shown 
in Fig. 19, use the Alexander & Ferguson (1994) and Iglesias & Rogers (1996) opacities, and 
the model atmospheres of Plez (1992) and Kurucz (1991). They treat convection with MLT and 
ctin = 1.6 to match the present-day Sun. 

These tracks behave in a similar fashion as the PS99 tracks. In particular, the position of the 
V1174 Ori primary is precisely matched by the theoretical 1.01 M Q track, but the agreement with 
the location of the V1174 Ori secondary is extremely poor. Indeed, the theoretical 0.73 M track is 
a full ~ 500 K too hot compared to the observed position of the secondary star. These tracks, like 
the others considered above, appear to be too compressed in T e g to simultaneously fit the primary 
and secondary of VI 174 Ori. 

As with the Palla & Stahler (1999) models, these tracks yield non-coeval ages for the two 
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components of V1174 Ori. Once again, the inferred ages for the primary and secondary stars differ 
by a factor ~ 2: in this case, ~ 7 Myr and ~ 3 Myr for the primary and secondary, respectively. 
But as with the PS99 models, despite non-coevality in the H-R diagram, these models produce 
isochrones in superb agreement with the observations in the M-R plane (7 Myr). 

6.4. MDKH03 

Finally, we consider the new PMS evolutionary tracks of Montalban et al. (2003), descended 
from the earlier models of D'Antona & Mazzitelli (1997). These authors have produced a variety 
of models that treat convection via MLT or FST, with the non-gray atmospheres from Hauschildt, 
Allard, & Baron (1999) or Heiter et al. (2002), respectively. For the models using the NextGen 
atmospheres of Hauschildt, Allard, & Baron (1999), the authors compute families of models with 
different values of T p h- Here we consider three models as listed in Table 11. These models are 
compared with the measured positions of the components of V1174 Ori in Figs. 20-21. 

The MLT models with NextGen atmospheres and aj n = 1.0 are most comparable to the 
BCAH98 tracks in Fig. 16. As the only substantial difference between these models is the depth of 
the outer atmosphere (r p h = 3 for MDKH03 tracks vs. 100 for BCAH98 tracks), these tracks yield 
very similar results. The only notable difference appears to be in the rate at which the models 
contract, with the r p h = 3 tracks yielding slightly younger ages, as seen in the M-R plane. 

In contrast, both the MLT and FST models, with the ATLAS9 atmospheres of Heiter et al. 
(2002), show excellent agreement with the position of the VI 174 Ori in the H-R diagram, with 
the FST models being ~ 100 K hotter. Unfortunately, these models do not permit a comparison 
with the V1174 Ori secondary as the Heiter et al. (2002) atmospheres do not extend below 4000 K. 
However, the general behavior of these tracks in the H-R diagram suggests that the 0.73 M track 
could intersect the position of the secondary if the calculation were extended to lower T c g . 

Indeed, these ATLAS9-based models differ qualitatively from all of the other models considered 
above in terms of the general shape of the tracks. These models exhibit a dramatic sweep to cooler 
T e ff, and then back to hotter T e g, as the stars descend their tracks, particularly for M < 0.8 M . 
The effect is to spread the mass tracks over a larger range of T e g at ages between 3-30 Myr, possibly 
consistent simultaneously with the primary and secondary of V1174 Ori. 

6.5. Other mass measurements of PMS stars 

As we have seen, the empirical mass measurements of V1174 Ori provide important constraints 
on the absolute mass calibration of PMS evolutionary tracks, both in the H-R diagram and in 
the M-R plane, between 0.7 M Q and 1.0 Mq. In this section, we combine our measurements of 
V1174 Ori with extant mass measurements of other PMS stars in order to extend our analysis to 



-31 - 



a larger range in stellar mass and to compare our findings among stars in both binary and single 
systems. 

Empirical mass determinations exist in the literature for twelve 15 other PMS stars — 8 in binary 
systems, 4 single — with masses accurate to better than ~ 10% and that are within the range covered 
by the PMS tracks considered above (< 1.6 M ). These are: TY CrA B with M = 1.64 ± 0.01 M 
(Casey et al. 1998); AK Sco A and B, both with M = 1.35±0.07 M (Alencar et al. 2003); EK Cep 
B with M = 1.12±0.01 M (Popper 1987); NTT 045251+3016 A and B with M A = 1.45±0.19 M 
and M B = 0.81+0.09 M (Steffen et al. 2001); RXJ0529.4+0041 A and B with M A = 1.25±0.05 M 
and M B = 0.91 ± 0.05 M (Covino et al. 2000); Lk Ca 15 with M = 0.97 ± 0.03 M (Simon, 
Dutrey, k Guilloteau 2000); GM Aur with M = 0.84 ± 0.05 M (Simon, Dutrey, & Guilloteau 
2000); DL Tau with M = 0.72 ± 0.11 M (Simon, Dutrey, & Guilloteau 2000); and DM Tau with 
M = 0.55 ± 0.03 M (Simon, Dutrey, & Guilloteau 2000). 

In Figs. 22-26 we show all of these stars together on H-R diagrams with the different PMS 
tracks discussed above. Binary and single stars are plotted separately in each figure. This ensemble 
of T e ff , luminosities, and the uncertainties therein are taken from the compilation of Luhman et 
al. (2003), which used the same temperature scale used in this study. To this compilation we have 
added the secondary of TY CrA (Casey et al. 1998), the two (identical) components of AK Sco 
(Alencar et al. 2003), and the two components of V1174 Ori. Thus the measurements represented 
in Figs. 22-26 are all based on a consistent T e g scale. 

We note that the listed uncertainties in the masses of the single stars do not include systematic 
uncertainties on the distances to those stars, which can amount to an additional ~ 15% uncertainty 
on the masses (Simon, Dutrey, & Guilloteau 2000). We note also that the meanings of the error 
bars in L and in T e g take on somewhat different meanings for the binary and single stars. For the 
single stars and the components of the astrometric binary NTT 04251+3016, these quantities are 
measured independently, whereas for stars in eclipsing binaries these quantities are coupled because 
L is derived from T e g and the stellar radii. Finally, we note that for the binary stars the individual 
masses are coupled to one another by virtue of the measured mass ratio; adjusting the mass of 
one component necessarily adjusts that of the other. These considerations constrain the extent to 
which individual objects can be "fine-tuned" into agreement with the models. 

We do not find a set of tracks that reproduces the observed locations of these stars simultane- 
ously at both high and low masses, and in both binary and single systems. 

Consider the BCAH98 tracks with ai n = 1.0 (Fig. 22), which have the coolest temperatures of 
all the tracks we have examined. These tracks agree with most of the binary stars, so long as the 
masses of individual objects are in some cases adjusted within their uncertainties. For example, 



15 We exclude from this discussion RS Cha A and B (M — 1.86 M Q and 1.82 M respectively), because of their 
high masses. We also exclude BP Tau (M = 1.32 Mq) because of the large uncertainty on the mass determination 
(±0.2 Mq) and the large uncertainty on its distance (Simon, Dutrey, & Guilloteau 2000). 
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for the tracks to be consistent with the primary of NTT 045251+3016, its mass must be adjusted 
downward by la. Of course, doing this requires that the mass of the secondary also be adjusted 
downward by la, implying a true mass of 0.72 Mq for the secondary. Even so, the tracks miss 
this star by ~ 1.5a. The greatest challenge to these tracks among the binaries is the 1.01 M Q 
primary of VI 174 Ori, which has the most accurate mass determination of the lot. This failure to 
match the VI 174 Ori primary is amplified by the position of LkCa 15; with a very similar mass of 
0.97 ± 0.03 M , it is located at almost precisely the same position in the H-R diagram. And like 
V1174 Ori A, LkCa 15 is not matched by these models. 

Indeed, none of the single stars are matched by these tracks, the tracks being systematically 
too cool. The single stars are better matched by the BCAH98 a in = 1.9 tracks (Fig. 23), which are 
~ 200 K warmer, a conclusion also reached by Luhman et al. (2003). However, these tracks have 
difficulty when compared against the binaries, especially for those with M < M . 

The other models show similarly mixed success in reproducing the observations. The PS99 
tracks (Fig. 24) perform comparably to the BCAH98 a m = 1.9 tracks, matching the positions of all 
of the single stars and, with the notable exception of NTT 045251+3016 A, matching the binaries 
with M > M (including the primary of V1174 Ori) very well. However, the lower-mass binary 
components — including the secondary of V1174 Ori — are poorly matched by these tracks. The 
results for the SDF00 models (Fig. 25) are very similar as well, though with the warmest and most 
closely spaced tracks of those considered here, the discrepancies among the lowest-mass stars are 
particularly severe, even for the single star DM Tau. Finally, a full evaluation of the MDKH03 
MLT ATLAS9 tracks (Fig. 26) is not possible because they do not extend below 4000 K. However, 
based on the qualitative behavior of these tracks it appears possible that they could perform well 
over a broad range of masses and for both binary and single stars. The lower temperatures of 
these tracks means that they agree with many of the higher-mass binary stars (though again with 
NTT 045251+3016 A being a clear outlier), and with their horizontal sweep in T e g — effectively 
stretching the T e g scale for stars with M < M — may be able to simultaneously match the positions 
of the lower-mass binaries and single stars. 

For completeness, we show in Figs. 27-28 the mass-radius relation for the three PMS binaries 
for which empirical mass measurements exist for both components (V1174 Ori, NTT 045251+3016, 
and RXJ0529.4+0041). With the relatively large error bars on the components of NTT 045251+3016 
and RXJ0529. 4+0041, we can draw limited conclusions here. All of the models show qualitative 
agreement with the general slope of the mass-radius relationship. However, the agreement does 
appear to be best among the tracks with low values of G^n; 

with the BCAH98 a m — 1.0 tracks and 
the PS99 a m = 1.5 tracks providing the best fit to the observations. The SDF00 a- m = 1.6 tracks 
perform reasonably well, whereas the tracks with large a- m (e.g. BCAH98 with a- m = 1.9) are least 
consistent with these measurements. 

It is valuable and instructive to consider what insights can be gained from the observations 
alone, independently of the models. First of all, it appears clear that the masses of the components 
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of NTT 045251+3016 are over-estimated. The primary star is at a very similar position in the H-R 
diagram to that of the 1.01 M primary of V1174 Ori, despite having a measured mass that is 
considerably larger. It does have a somewhat higher luminosity, so it is possible that it does indeed 
have a slightly higher mass than does V1174 Ori A. This conclusion is bolstered by the secondary 
of NTT 045251+3016; its position in the H-R diagram, suggesting a mass slightly less than the 
0.73 M secondary of V1174 Ori, belies its measured mass of 0.81 M Q . That NTT 045251+3016 
may be an outlier may not be surprising given the large uncertainties inherent to the analysis of 
this astrometric binary, for which only ~ 1/4 of the orbit has been mapped (Steffen et al. 2001). 
Indeed, simply adjusting the mass of the NTT 045251+3016 primary down by 1-2<j, while keeping 
the mass ratio fixed, would appear to bring the measurements of NTT 045251+3016 and V1174 Ori 
into closer empirical agreement. 

Similar concerns arise with the components of the eclipsing binary RXJ0529. 4+0041. Here, the 
position of the primary star in the H-R diagram coincides closely with that of TY CrA B, despite 
the latter having a precisely determined mass that is 30% larger. At the same time, however, 
the secondary of RXJ0529. 4+0041 is located beween the two components of V1174 Ori, which is 
appropriate for its measured mass intermediate beween those of VI 174 Ori A and B. This suggests 
that the primary (and perhaps the secondary as well) are simply over-luminous, which would tend to 
implicate the stellar radii in this case since, again, the T e g- of the secondary appears approximately 
correct relative to the components of V1174 Ori. Upon re-examination of the Covino et al. (2000) 
results for RXJ0529. 4+0041, it appears that the uncertainties they report for the luminosities 
are underestimated. For example, the fractional uncertainty reported by them for the luminosity 
of RXJ0529. 4+0041 A (8%) is smaller than the fractional uncertainty they report (11%) for the 
radius alone. Including the uncertainty in T e fj- as well gives a more appropriate uncertainty in the 
luminosity of ~ 15%. 

In some cases, the observations place objects of similar masses at similar positions in the H-R 
diagram, bolstering their credibility. For example, the primary of V1174 Ori and LkCa 15, both with 
M « 1.0 M , have very nearly the same T e g and L. The same is true for RXJ0529. 4+0041 B and 
GM Aur, both with M « 0.87 M . Even though these latter measurements have somewhat larger 
uncertainties, the lower luminosity of RXJ0529. 4+0041 B relative to (the presumably younger) 
GM Aur in effect provides an empirical mass track for M « 0.87 M between ~3-10 Myr. 

But there are also cases where the data seem at odds with one another. For example, V1174 Ori 
B and DM Tau occupy similar positions in the H-R diagram despite having masses of 0.73 ± 
0.01 M and 0.55 ± 0.03 M Q , respectively. This is worrisome. One way to reconcile this apparent 
contradiction is to have mass tracks, like the MLT ATLAS9 tracks of Montalban et al. (2003), that 
are double-valued in T e g at low masses, allowing stars with different masses to have similar T c g 
at certain times. While these tracks do not presently extend to the masses of VI 174 Ori B and 
DM Tau, inspecting these tracks in Fig. 26 one can imagine that these stars' positions would be 
well matched by their respective tracks, particularly with DM Tau being somewhat more luminous 
than VI 174 Ori B. 
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Another possible explanation for the discordance between these measurements is that the 
uncertainties in the masses for the single stars are underestimated. Indeed, in contrast to the 
distance-independent masses derived for V1174 Ori, the masses determined for the single stars (from 
spatially resolved disk rotation profiles; Simon, Dutrey, & Guilloteau (2000)) depend directly upon 
the distance assumed. As noted by Simon, Dutrey, Sz Guilloteau (2000), the distance uncertainty 
to any particular star in Taurus may be as large as ~ 15%, manifesting itself as a systematic error 
in the dynamical mass determination. 

DL Tau may serve as a case in point. Its measured mass of 0.72 M Q is similar to the mass 
of 0.73 M we have measured for the secondary of V1174 Ori. Yet DL Tau is 450 K hotter than 
V1174 Ori B, suggesting a mass more comparable to that of V1174 Ori A. An important additional 
clue to the relative masses of these two objects is provided by their relative Li abundances. In 
Fig. 29 we show the Li abundances, log ./V (Li), for all of the PMS stars discussed above for which 
Li measurements exist in the literature. These stars exhibit a clear correlation of log iV(Li) with 
stellar mass: Stars with M > 1.2 M show undepleted levels of Li, while less massive stars show 
a trend of monotonically decreasing Li abundance with decreasing stellar mass. DL Tau is a clear 
outlier in this trend. Its measured logiV(Li), even with its large uncertainties, is too high for its 
mass. In contrast, V1174 Ori B follows the empirical trend nicely. The observed Li abundance 
of DL Tau implies a lower limit to its mass of > 0.9 Mq, consistent with its location in the H-R 
diagram relative to V1174 Ori A. Whether the masses of V1174 Ori B and DM Tau can be similarly 
reconciled is unknown; there is no published Li abundance for DM Tau. 

The observed pattern of Li depletion with stellar mass provides more general constraints on 
PMS evolutionary tracks, particularly with respect to the efficiency of convection. As discussed 
in detail by D'Antona & Montalban (2003), models with more efficient convection (i.e. larger a) 
result in dramatically more Li depletion than do models with less efficient convection. 

Consider a star like V1174 Ori A, having M = 1.0 M Q and situated near the bottom of its 
Hayashi phase of evolution. Assuming an initial Li abundance of logiV(Li)= 3.3, a model with 
cti n = 1.9 will have logiV(Li) ~ 2.5 whereas a model with a more modest a- m = 1.0 will have 
log N(Li) ~ 3.2 (D'Antona & Montalban 2003). The primary of V1174 Ori is more consistent with 
the low Li depletion predicted by the a- m = 1.0 model. At lower masses, the depletion predicted 
for models with highly efficient convection becomes ever more severe. A star with M = 0.7 M Q 
is predicted to have logiV(Li) ss —1.5 by the a- m = 1.9 models, completely inconsistent with the 
modest depletion levels seen in Fig. 29, which agree with the predicted levels of depletion for 
a m = 1.0 models (D'Antona k, Montalban 2003). Thus the Li depletion measurements strongly 
favor models with low a; n . 

At low masses, the Li depletion levels predicted by the models also become strongly dependent 
upon the parameter r p h. Larger values of r p h result in a larger fraction of the over-adiabatic 
region being computed with the low value of a adopted in the model atmosphere. This produces 
a reduced temperature at the base of the convection zone, resulting in lower levels of Li depletion. 
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A more detailed analysis of the observed Li abundances as a function of stellar mass is beyond 
the scope of this paper. But coupled with the precise mass determinations for PMS stars, such 
as we have determined here, the clear pattern of Li depletion with stellar mass should in principle 
provide important additional constraints on various key parameters of PMS evolutionary models. 
Moreover, as we have seen here, Li abundances can provide valuable, independent mass constraints 
as well when systematic uncertainties (e.g. distance) may limit the accuracy of the dynamical 
measurements. 

7. Summary and Conclusions 

We have identified the star V1174 Ori as a previously unknown, double-lined, spectroscopic, 
pre-main-sequence (PMS) eclipsing binary system in the Orion star-forming region. 

From high-cadence, multi-epoch, multi-band, time-series photometry spanning 2900 days, we 
find the period of the system to be 2.634727 days. From multi-epoch, high-resolution spectroscopy, 
we determine a precise, double-lined orbit solution which constrains the absolute masses of the 
primary and secondary components to be 1.01 M and 0.73 M , respectively, accurate to ~ 1%. 
Synthetic light-curve modeling and analysis provides absolute stellar dimensions. We find the stars 
to have radii of 1.34 R Q and 1.06 R , respectively, again with uncertainties of ~ 1%. These 
empirical mass and radius measurements are distance-independent. The secondary is thus the 
lowest-mass PMS star yet discovered in an eclipsing binary system, and the first sub-solar PMS 
star with an empirically determined mass accurate to better than a few percent. 

We determine the distance to the system to be ~ 400 pc which, together with the system's 
center-of-mass radial velocity, makes V1174 Ori a likely member of the Ori OBlc region. It is 
therefore likely slightly in the foreground of, and slightly older than, the Orion Nebula Cluster 
(ONC) proper. Given the spread of available estimates for the bolometric magnitude of the Sun, 
we cannot currently constrain the absolute distance to VI 174 Ori more stringently. 

Analysis of the system colors suggests the presence of an unseen tertiary companion in a wide 
orbit. Follow-up photometry in the near-IR is underway, and should reveal substantial third light, 
which we predict to be ~ 30% in J. 

We compare the empirically determined stellar dimensions with a variety of PMS stellar evo- 
lution models in both the H-R diagram (T c g, L plane) and the mass-radius plane. In all cases, 
the models yield an age for V1174 Ori of ~ 5-10 Myr, consistent with membership in the slightly 
evolved Ori OBlc region, just in the foreground of the ONC. 

While all of the models predict the slope of the mass-radius relationship reasonably well, those 
with less efficient convection (i.e. small convection parameter a) are the most consistent with our 
measurements, even though larger values of a are required to match the present-day Sun. While 
some models are consistent with the mass of the V1174 Ori primary, and others are consistent 
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with the mass of the V1174 Ori secondary, none of the models correctly predict the masses of both 
stars simultaneously. This apears to be the result of the theoretical mass tracks between 0.7-1.0 
M Q being too compressed in T e g relative to the components of V1174 Ori. Recent models by 
Montalban et al. (2003) employing the ATLAS9 atmospheres of Heiter et al. (2002) seem to show 
a larger spread of T e g in this mass range, and could possibly fit both components of VI 174 Ori, 
but these models do not presently extend to masses low enough to test against the V1174 Ori 
secondary. 

Potential problems with the T e g scale are also implicated by a curious effect in which the 
theoretical tracks and the empirical measurements translate differently between the mass-radius 
plane and the H-R diagram. Some of the tracks show the components of V1174 Ori to be coeval in 
the H-R diagram but non-coeval in the mass-radius plane, and vice-versa for others. 

Analysis of the H-R diagram positions of all PMS stars with empirical mass measurements 
yields mixed results for all of the PMS models we examined, with no one set of models clearly 
superior to the others. However, the ensemble of observations appear to be least compatible with 
the tracks of Siess, Dufour, & Forestini (2000), which have the hottest and most compressed T e g 
scale of any of the tracks we examined. Again, the ATLAS9-based models of Montalban et al. (2003) 
show promise of being able to reproduce most of the empirical measurements, and we encourage 
the extension of these models to lower T e g. 

We have carefully scrutinized the ensemble of empirical PMS mass measurements for consis- 
tency with one another, independent of the models. Encouragingly, we find instances of superb 
agreement between independent mass determinations, such as V1174 Ori A and LkCa 15 (a single 
star in Taurus), which have identical masses and are at identical positions in the H-R diagram. 
But, disconcertingly, we also find instances of apparent contradiction, such as V1174 Ori B and 
DM Tau, which are at similar (though not identical) positions in the H-R diagram yet have very 
different masses. 

These apparently discordant empirical measurements may be reconcilable, however. PMS 
evolutionary models that are double- valued in T e g , such as the ATLAS9-based models of Montalban 
et al. (2003), may in fact predict the H-R diagram positions observed for these stars. In some 
cases, such as the astrometric binary NTT 045251+3016 and the eclipsing binary RXJ0529. 4+0041, 
we reconcile the disagreements by adjusting the masses of the stars within the uncertainties. In 
another case — the single star DL Tau — we find compelling evidence in the form of an anomalous 
Li abundance that the reported mass is incorrect, most likely the result of a systematic error in the 
assumed distance. 

Indeed, we find that the ensemble of PMS stars with empirical mass determinations possess 
a remarkably well-defined pattern of Li depletion as a function of stellar mass. The mild amount 
of depletion observed in the lowest-mass star, V1174 Ori B, strongly constrains the efficiency of 
convection that can be acting in the interiors of these stars, and argues strongly in favor of models 
with small convection parameters, a. Thus, as argued by D'Antona & Montalban (2003), Li 
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abundance determinations in PMS stars promise to provide an important observational window 
into the interior structure of PMS stars. And, as we have shown in the case of DL Tau, Li can 
provide a valuable means for triangulating against the dynamical masses of these stars in the 
presence of systematic errors. 
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Table 7. Radial velocities of HD 26162 



Radial Velocity 

(km/s) 



1 


25.03 


± 


0.08 


2 


25.33 


± 


0.08 


3 


25.17 


± 


0.07 


4 


24.76 


± 


0.07 


5 


23.61 


± 


0.08 


6 


25.20 


± 


0.07 


7 


25.35 


± 


0.07 


9 


25.51 


± 


0.07 


10 


25.38 


± 


0.07 


11 


24.96 


± 


0.08 


12 


24.74 


± 


0.09 


13 


24.33 


± 


0.07 


14 


25.24 


± 


0.07 


15 


25.49 


± 


0.08 



a Scc Tabic 3. 
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Table 8. Final light curve solution for VI 174 Ori 



parameter 


value 


parameter 


value 




86.974 ± 0.035 


ta pole 


0.1384 ± 0.0006 


n A 


7.942 ± 0.026 


point 


0.1392 ± 0.0006 


n B 


7.705 ± 0.022 


side 


0.1387 ± 0.0006 


TA (K) 


4470. (fixed) 


back 


0.1391 ± 0.0006 


TB (K) 


3615.0 ± 9.5 


re pole 


0.1101 ± 0.0012 


is(J) (%) 


9.90 ± 0.71 


point 


0.1106 ± 0.0012 


L 3 (V) (%) 


0.00 (fixed) 


side 


0.1103 ± 0.0012 


^3 (2/) (%) 


1.73 ± 0.65 


back 


0.1106 ± 0.0012 


L 3 (6) (%) 


0.03 ± 0.65 


Fout,A 


0.4247 ± 0.0012 


L 3 («) (%) 


2.47 ± 0.66 


Fout,B 


0.4548 ± 0.0014 


fratio.A 


1.017 ± 0.026 






fratio.B 


1.001 ± 0.079 







Note. — WD solutions for V1174 Ori, based on the observa- 
tions from Jan 2001. Parameters like the limb darkening coef- 
ficients, gravity darkening exponents, and surface gravities were 
kept consistent during iterations. The solutions were attained 
by combining both the least-squares and the SIMPLEX meth- 
ods. The errors quoted are the internal least-squares estimations, 
which are typically around 3-5 times larger then those estimated 
from the vertices of the SIMPLEX hypcrpolyhcdron after conver- 
gence. The third light contributions are given as percentages of 
the eclipsing system total light at the normalization phase (0.25). 
The radii resulting from the gravitational effective equipotcntials 
(in terms of the orbital separation) and the fill-out factors, F ou t 
(Mochnacki 1984), are given, which indicate that the eclipsing 
components are well inside their Roche lobes. The errors quoted 
on the radii and fill-out factors represent the deviations caused 
by the uncertainties in the pscudopotcntials and in the mass ra- 
tio. The rotation rates, calculated from the solutions and the 
measurements of §4.3.1 are shown, justifying the assumption of 
synchronism with the orbit. 
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Table 9. Spot parameters used in light curve solutions 



epoch 


bands 


spot /star 


co-latitude (°) 


longitude (°) 


radius (°) 


T f 


actor 


Dcc/1994 


I 


1 


/ 


A 


69. 


± 


13. 


140.5 ± 8.6 


7.2 ± 1.6 


0.69 


± 0.24 






1 


/ 


B 


70.4 


± 


7.3 


25.3 ± 5.5 


11. ± 11. 


1.44 


± 0.32 


2MASS 


JHK 


1 


/ 


A 


21. 


± 


16. 


101.9 ± 9.6 


27. ± 5.3 


0.88 


± 0.18 






1 


/ 


B 


23. 


± 


13. 


336. ± 21. 


48. ± 4.2 


0.61 


± 0.14 


Jan/2001 


vbyVI 


1 


/ 


A 


18.2 


± 


7.4 


78. ± 10. 


15.2 ± 4.6 


1.58 


± 0.74 






2 


/ 


A 


16.2 


± 


5.7 


110.8 ± 6.2 


26.0 ± 4.7 


0.76 


± 0.25 






3 


/ 


A 


13.34 


± 


0.69 


159.6 ± 1.3 


31.2 ± 1.5 


0.508 


± 0.042 






1 


/ 


B 


155. 


± 


10. 


97.1 ± 5.5 


30.2 ± 5.1 


0.76 


± 0.12 






2 


/ 


B 


67. 


± 


19. 


246. ± 10. 


17.4 ± 4.6 


1.18 


± 0.11 






3 


/ 


B 


76. 


± 


17. 


297.4 ± 5.6 


26.0 ± 5.6 


0.662 


± 0.091 


Nov/2001 


I 


1 


/ 


A 


48. 


± 


10. 


136.8 ± 9.6 


21.4 ± 5.2 


0.74 


± 0.41 






1 


/ 


B 


70. 


± 


17. 


134.6 ± 9.7 


12.4 ± 6.6 


1.45 


± 0.91 


Nov/2002 


BVI 


1 


/ 


A 


22.41 


± 


0.42 


93.3 ± 7.2 


31.6 ± 1.2 


0.851 


± 0.057 






1 


/ 


B 


16.35 


± 


4.1 


306.1 ± 9.1 


56.1 ± 5.7 


0.69 


± 0.11 






2 


/ 


B 


21.99 


± 


4.2 


348.2 ± 9.8 


44.7 ± 8.4 


0.605 


± 0.095 



Note. — The characteristics of the star spots used in the solution of the light curves from Table 1 
and for the JHK light curves from 2MASS (one point obtained at March 19th, 1998 and 23 points 
from February 6th to March 28th, 2000). The spots are listed for each component in each epoch 
with increasing longitudes. As in Table 8, the errors quoted are those internal to the least squares 
method. The standard deviation of the corresponding variables of the SIMPLEX method, after 
convergence, are much smaller: of the order of 10~ 3 degree for the angular parameters and of 10~ 5 
for the temperature factor for both components. Tf actor gives the spot temperature in units of T e ff . 



Table 10. Fundamental stellar properties of V1174 Ori 





Primary 


Secondary 


M (M Q ) 


1.009 ±0.015 


0.731 ±0.008 


R (R Q ) 


1.339 ±0.015 


1.065 ±0.011 


T efI (K) 


4470^120 


3615^ 100 c 


logL (L Q ) 


-0.193 ±0.048 - 


-0.761 ±0.058 


logs 


4.19 ±0.01 


4.25 ±0.01 


v Iot (km/s) 


26.2 ±0.6 


20.5 ±1.6 


A v 


0.32 




Dist. (pc) 


419 ±21 



a From measured spectral type of K4.5 ±0.5, and 
SpT-T off conversion of Schmidt-Kalcr (1982). 

b Implies a spectral type of Ml. 5. 

c Note that the internal uncertainty on T c g for the 
secondary is minuscule. The uncertainty we assign 
here is inherited from the uncertainty in the deter- 
mination of T e g for the primary, keeping the T e g 
ratio fixed at the value determined in §4.4. 
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Table 11. Parameters of PMS Evolution Models 



Authors 


Atmosphere 


Convection 








Matches Sun? 


MDKH03 


ATLAS9 


FST 






10 


Y 




ATLAS9 


MLT 


0.5 


2.3 


10 


Y 




Next Gen 


MLT 


1.0 


1.0 


3 a 


N 


BCAH98 


Next Gen 


MLT 


1.0 


1.0 


100 


N 




Next Gen 


MLT 


1.0 


1.9 


100 


Y 


PS99 


gray approx. 


MLT 




1.5 


2 

3 


N 


SDFOO 


Plez, Eriksson 


MLT 




1.6 


10 


Y 



a Thesc authors have also computed NcxtGcn-based models with T p h = 100, but 
these are equivalent to the BCAH98 tracks. 
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Fig. 1. — Keck HIRES spectrum of V1174 Ori, in the vicinity of the Li A6708 line, is shown as 
a solid curve. The dashed curve is a spectrum of the V1174 Ori primary, obtained by averaging 
together all of the HET HRS spectra. The Li line from the secondary (labeled Li B) is clearly 
present, in this case blueward of the primary's Li line (labeled Li A) and partially blended with an 
Fe line from the primary. The emission feature near A6716 is nebular [S II]. 
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Fig. 2. — The vbyVI light curves observed in Jan 2001 together with the theoretical light curves 
synthesized with a Wilson-Devinney based algorithm (see Tables 8 and 9). For clarity, the light 
curves were shifted by the amount indicated in the figure, but the relative scale is the same for all 
passbands. Note the changes in the depths of both eclipses with wavelength. 
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Fig. 3. — The O — C residuals for the solutions of Fig. 2 and Tables 8 and 9. The vertical bar at 
the left of each panel is the standard deviation of the residuals. 
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Fig. 4. — The IV B light curves observed in 2002 and the / light curve obtained in 2001 (Table 1) 
together with the theoretical light curves obtained with the WD algorithm (Tables 8 and 9). The 
light curves are shifted by the amount indicated in the figure for clarity, but the relative scale is 
the same for all passbands. Note the changes in the depth of both eclipses with wavelength and 
the clear change in the form of the light curve outside eclipse from 2001 to 2002. 
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Fig. 5. — The O — C residuals for the solutions of Fig. 4 and Tables 8 and 9. The vertical bar at 
the left of each panel is the standard deviation of the residuals. 
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Fig. 6. — The / light curve obtained in Nov 1994 (Table 1) and the JHK light curves obtained by 
Carpenter, Hillenbrand, & Skrutskie (2001) together with the theoretical light curves synthesized 
by the WD algorithm (Tables 8 and 9). The light curves are shifted by the amount indicated in 
the figure for clarity, but the relative scale is the same for all passbands. Note the changes in the 
depth of both eclipses with wavelength and the clear change in the form of the light curves outside 
eclipse from 1994 to 2000. 
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Fig. 7. — The O — C residuals for the solutions of Fig. 6 and Tables 8 and 9. The vertical bar at 
the left of each panel is the standard deviation of the residuals. 
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Fig. 8. — The O — C phase residuals of the eclipse timings of V1174 Ori are shown based on a 
best-fit period of 2.634727 days. Filled circles represent primary eclipses, open circles represent 
secondary eclipses. See Table 4. 
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Fig. 9. — Example of our procedure for measuring radial velocities from cross-correlation, (a) Shown is the cross- 
correlation function from a single spectral order in the spectrum of UT 2002 Jan 28, in which we detect both the 
primary and the secondary peak. In a typical spectrum, most of the orders yield a velocity measurement for the 
primary, whereas only a few orders produce an unambiguous secondary peak, such as the one seen in this example, 
(b) Shown as a function of correlation peak height are the radial velocities from the 24 orders in the spectrum of UT 
2001 Dec 18 that give radial velocities (km/s) for the primary star. We use only velocities with correlation heights 
greater than 0.6 (asterisks), which we found empirically yield more reliable velocities. The horizontal lines indicate 
the mean velocity of the primary (72.91 km/s) and the standard deviation of the mean (0.19 km/s); see Table 5. 
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Fig. 10. — Illustration of distortions in observed radial velocities due to stellar occultation and the 
presence of surface spots. Distortions in the primary (secondary) star's radial velocities are shown 
by the solid (dashed) line. Note the strong distortion effects due to occultation at ±0.05 phase 
around each eclipse (±3.7 km/s and ±4.8 km/s for the primary and secondary, respectively). 
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Fig. 11. — Orbit solution for VI 174 Ori. Radial velocity measurements (see Table 5) for the primary 
(secondary) star are shown as filled (open) points. The solutions for the primary and secondary 
stars are shown as solid and dashed lines, respectively. 
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Fig. 12. — Our assessment of mechanical stability of the HET HRS. (a) Shown as a function of heliocentric Julian 
date are the radial velocities of the radial- velocity standard HD 26162 (see Table 5). Circled points represent spectra 
from the same nights that are used in our orbit solution of V1174 Ori. The horizontal dashed line indicates the 
mean of these points, which show a small systematic offset of 0.36 km/s, likely due to spectral mismatch between 
HD 26162 and the standard star (HD 18884) used in our orbit solution of V1174 Ori. (b) These same radial velocities 
of HD 26162 are shown as a histogram (open), with the circled points in (a) represented by the filled histogram. 
These latter measurements show an r.m.s. scatter of 0.26 km/s. 
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Fig. 13. — Shown is the UT 2002 Feb 07 spectrum of V1174 Ori (second from top) compared 
with three spectral standards of type K3 (bottom), K4 (second from bottom), and K5 (top). The 
spectral region shown includes the temperature-sensitive lines of Ni I at A6108 and V I at A6112, 
one of four line pairs that we use to determine the spectral type of the primary of V1174 Ori. 
Note that in this spectrum the primary star's V I is slightly contaminated by the secondary star's 
Ni I line. Also visible here are the temperature-insensitive lines of Ca II at AA6102,6122. In this 
spectrum of VI 174 Ori, the secondary star's A6102 Ca II line is visible as the weak absorption 
feature between the primary star's A6102 Ca II and A6108 Ni I lines. The A6102 Ca II lines in this 
spectrum indicate a primary-to-secondary flux ratio of ~ 6:1 for VI 174 Ori (see Fig. 14). 
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Fig. 14. — Determination of the primary-to-secondary flux ratio for V1174 Ori. The UT 2002 Feb 
07 spectrum (solid line) is compared with a synthetic spectrum produced by combining the K3 and 
K7 standards in a 6:1 ratio (dashed line). The standards were first rotationally broadened (see 
§4.3.1) and shifted to the appropriate radial velocities (see Table 5). 
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Fig. 15. — Determination of Ha equivalent widths for the two components of VI 174 Ori. Shown is 
the UT 2002 Feb 23 spectrum, plotted vs. heliocentric radial velocity. The strong emission line is 
nebular. Emission from the primary and secondary stars are visible just to the blue and red of the 
nebular line, respectively. 
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Fig. 16. — BCAH98 tracks with lmi*/H p = 1.0. (a) H-R diagram (logL vs. T cff ), with isochrones at 1, 3, 10, 30, 
and 100 Myr, and mass tracks for 0.6, 0.7, 0.8, 0.9, 1.0, and 1.1 M©. Mass tracks at 0.7 and 1.0 M©, corresponding 
approximately to the masses of V1174 Ori, are thicker than the others. Points mark the best-fit values of T e g 
and log L, with nested boxes indicating confidence regions about those values. The inner box about the primary 
corresponds to 0.5 spectral sub-type uncertainty in T e s (see §4.3.3) and la uncertainty on the stellar radius. The 
outer box corresponds to uncertainties twice as large in each of T e g and R, and should be taken as the region of high 
confidence for the primary. The nested boxes about the secondary are similar, but the range of T e g here is determined 
by keeping the secondary-to- primary T e a ratio fixed at the value determined from the light curve analysis (see §4.4). 
Note that the positions of the primary and secondary in this figure are highly correlated; adjusting the position of 
the primary necessarily adjusts that of the secondary, (b) Radius vs. mass in solar units. Points with error bars are 
the measurements for the two components of V1174 Ori; isochrones at 3, 10, 30, and 100 Myr are shown as dashed 
lines. 
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Fig. 18. — PS99 tracks. Symbols and lines are as in Fig. 16, except that (a) shows mass tracks at 
0.6, 0.7, 0.8, 1.0, and 1.2 Mq and isochrones at 1, 3, 10, and 30 Myr. Note that the mass track at 
0.7 M has been visually interpolated between the positions of the 0.6 and 0.8 M tracks. In (b) 
isochrones are shown at 3, 10, and 30 Myr. 
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Fig. 20. — Montalban et al. (2003) tracks with MLT convection (ai n = 1.0, r p h = 3) and NextGen 
atmospheres from Hauschildt, Allard, & Baron (1999). Symbols and lines are as in Fig. 16. 
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Fig. 21.— Montalban et al. (2003) tracks with (a) MLT convection (a in = 2.3) and (b) FST 
convection. Both use the ATLAS9 atmospheres from Heiter et al. (2002). Symbols and lines are as 
in Fig. 16. Note that the tracks are truncated below 4000 K as the Heiter et al. (2002) atmospheres 
do not extend to cooler temperatures. 
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Fig. 22. — Comparison of all PMS stars with accurate empirical mass determinations and M < 1.6 
M Q with the tracks of BCAH98, a- m = 1.0. (a) Stars in binary systems. Masses for each of these 
are indicated at the top of the figure. Mass tracks shown are for 0.6, 0.7, 0.8, 0.9, 1.0, 1.1, 1.2, and 
1.4 Mq. (b) Same as above, but for single stars. Note that the uncertainties on the masses of the 
single stars do not include distance uncertainties, which can amount to ~ 15% (Simon, Dutrey, & 
Guilloteau 2000). 
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Fig. 24. — Same as Fig. 22 but for PS99 tracks. Mass tracks shown are for 0.6, 0.7, 0.8, 1.0, and 
1.2 M m . 
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Fig. 27. — Mass-radius relationships (solar units) compared to BCAH98 tracks with ct\ n = 1.0 (top) 
and ai n = 1.9 (bottom). 




Fig. 28.— 



Same as Fig. 27 but for PS99 and SDFOO tracks. 
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Fig. 29. — Li abundances of PMS stars with empirical mass determinations, (a) Li abundance plotted as a function 
of stellar T e g (K). The two components of V1174 Ori are shown as filled circles, other stars are shown as open 
circles, (b) Li abundance plotted as function of empirical stellar mass (solar units). In both plots, the "cosmic" Li 
abundance of 3.31 is indicated by the horizontal line. Abundances are calculated from the observed Li equivalent 
widths, using the NLTE curve-of-growth calculations of Pavlenko & Magazzu (1996). References for the equivalent 
width measurements are as follows: VI 174 Ori (this study); RXJ0529.4+0041 (Covino et al. 2000); NTT 045251+3016 
A (Walter et al. 1988); TY CrA B (Casey et al. 1998); LkCa 15 (Martin et al. 1994); GM Aur (Basri, Martin, & 
Bertout 1991); and DL Tau (Basri, Martin, & Bertout 1991). 



